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© der „Zeitschrift für angewandte Mathematik und Mechanik“. 
ta zweijähriger Unterbrechung erscheint mit vorliegendem Heft 1 des Bandes 25/27 die „Zeit- 55 
schrift für angewandte Mathematik und Mechanik‘ wieder in der alten bewährten Anordnung 
"und Form. Das letzte Heft vor dem Zusammenbruch war das Doppelheft 5/6 des Bandes 24, das Ludwig 
Prandtl zum 70. Geburtstag gewidmet war. Heft 3 und 4 dieses Bandes waren zwar vorbereitet, konnten 
aber nicht mehr gedruckt werden. Dieser Band bleibt also unvollständig. er 
Künftig wird die Zeitschrift monatlich in Heften von 32 Seiten Umfang erscheinen. 
In den neuen Bänden will die Zeitschrift nicht nur, wie es in den letzten Jahrgängen der Fall 
war, Aufsätze in der Hauptsache aus den Gebieten der praktischen Analysis und der Mechanik fester 
und deformierbarer Körper bringen, sondern es sollen auch Arbeiten aus anderen Gebieten, z. B. der 
Elektrotechnik, der Thermodynamik, der Optik, der Geodäsie, der Statistik usw. aufgenommen werden, \ 
Jalls sie mathematisch interessante Gedankengänge enthalten. Die Zeitschrift soll wieder so umfassend 
werden, wie es die ersten Jahrgänge waren. In den letzien Jahren lag das Schwergewicht auf der 
Mechanik. In Zukunft sollen die Anwendungen der Mathematik m Vordergrund stehen. 
Der neue Verlag hat die Lizenz für die „Zeitschrift für angewandte Mathematik und Mechanik“ 
am 27.2.1947 erhalten. Seine umfangreiche wissenschaftliche Verlagsproduktion, vor allem auch 
auf verwandten Gebieten, gibt der „Zeitschrift für angewandte Mathematik und Mechanik“ das sichere 
Fundament für eine stete Aufwärtsentwicklung. gi 
Wir bitten num alle Fachgenossen des In- und Auslandes um ihre Mitarbeit, damit der Zeit- 
schrift im Interesse der internationalen Wissenschaft ihre alte-geachtete Stellung erhalten. bleibt. 
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HAUPTAUFSÄTZE 


Über die Anfangswertaufgabe in der Theorie 
der Rayleighwellen. 
Von Erwin Hardtwig in Stuttgart. 


Die Bewegungen, die von den Massenpunkten eines, den Halbraum erfüllenden elastischen (oder fast- 

- elastischen) Mediums ausgeführt werden, hängen davon ab, aus welchen Anfangsbedingungen heraus sie 

erfolgen. Die Bestimmung dieser Bewegungen wird für den Fall durchgeführt, daß Anfangsverschiebung 

und Anfangsgeschwindigkeit der Oberflächenpunkte vorgegeben sind. Der zwei- und dreidimensionale Fall 
werden getrennt behandelt. 

Les mouvements descrits par les points de masse d’un milieu elastique (ou presqu’ Elastique) remplissant 
le demi-espace dependent des conditions initiales d’apres lesquelles ils se produisent. La determination de 
ces mouvements est r&alisee pour le cas, ou deplacement initial et vitesse initiale des points superficiels sont 
indiques. Le cas de. deux et trois dimensions sont traites separement. 

T'he movements performed by the points of mass of an elastic (or nearly elastic) medium filling the half- 
room depend on the initial conditions from which they arise. T’he determination of these movements is per- 
formed for the case that initial displacement and initial velocity of the points of superficies are indicated. 
The two dimensional and three dimensional cases are treated. separately. 

Bu:keHuus, BEITONHAEMEIe MACCOBBIMU HYHKTAMH HATIONIHAIOIIETO 06, eM UMJIHHAP& 9NACTHIHOTO 
(HH HOYTH YTO 9NACTHIHOTO) MeHyMa, BABHCHT OT TOTO, 3 KAKHX TIEPBOHAMAJIBHBIX YCNOBHU 
OHH ucxonart. Onpexenernue aTHx ABHKeHUK IPOBONHTCH H& TOT CIIYYaH, UTO YKASAHLI HAUYAJIB- 
HOeE TepenBH:KeHHe U HAYAJIBHAA CKOPOCTB IYHKTOB IIOBEPXHOCTH. Cayyau ABOKHOTO H TPOAHOTO 
H3MePeHHA NOJIKHEI PACCMATPHBATECH B OTAENIBHOCTH. 


Von den drei Wellenarten, die in einem elastischen Körper bei entsprechender Anregung 
entstehen können, den longitudinalen Raumwellen, den transversalen Raumwellen und den 
Oberflächenwellen (Rayleighwellen) bieten zweifellos die letzteren die interessantesten Frage- 
stellungen sowohl praktischer wie theoretischer Art. Abgesehen davon, daß sie nicht in jedem 
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beliebig geformten Körper auftreten können, legen sie auch der Wahl der Anfangsbedingungen, 
aus denen heraus das System in freie Schwingungen übergeht, gewisse Einschränkungen auf. 

Die folgenden Zeilen befassen sich mit einer dahinzielenden Fragestellung. Um sie zu formu- _ 
lieren, bedienen wir uns, wie üblich, eines rechtwinkligen, kartesischen Koordinatensystems 
x, y, 2, das wir so legen, daß der unter der &-y-Ebene liegende Halbraym vom Medium erfüllt 
ist. Die positive z-Achse weise nach unten, also gegen das Medium hin, von dem wir annehmen, 
daß es die sich nur wenig zeitlich und räumlich ändernde Dichte o besitze. Die drei längs den 
Koordinatenachsen genommenen Komponenten des Verschiebungsvektors 3 seien u, v, w. 
Für den Spannungstensor mit den Komponenten 9,9% Pay - - machen wir, um möglichst die 
Allgemeinheit zu wahren, einen Ansatz, der sowohl den elastischen als auch den viskosen Eigen- 
schaften des Mediums gerecht werden soll und der in der Literatur wiederholt dort ange- 
wendet wurde, wo es sich darum handelte, das physikalische Verhalten von nicht vollkommen 
elastischen Körpern zu beschreiben). Wir verstehen unter A, u die Lameschen Konstanten | 
und unter A’, « zwei weitere konstante, Zahlen, die die viskosen Eigenschaften des Mediums 
charakterisieren sollen. Der Ansatz soll dann ‚so lauten: 


| EL ‚ou ER Ei 
Du =A0 Hau" +2 RO +2u au, Pyy Tree: Pız =: ” 
ov ou 0 [dv ou ’ 2 E 
el ea rn a Ehe ie a EN) 
Pay (+55) m +) PL Bus ” 


Der so angenommene Spannungstensor steht mit den zeitlichen Veränderungen des Ver- 
schiebungsvektors in Beziehungen, die, vektoriell geschrieben, sich durch die Differential- 
gleichung 

025 


0 = = 
U +(4 + u) grad div 3 +3, [#48 +(A + u) grad div 5] a 


für 3 ausdrücken lassen. Dazu kommt noch, als Oberflächenbedingung, die Forderung, daß 
die zur Oberfläche z=0 des Mediums normalen Spannungskomponenten für 2=0 gleich Null 
werden: 


BR ee 


Aus (I) ist ersichtlich, daß man die Gleichungen. der Elatizitätstheorie erhält, wenn man 
A=u =0 setzt und jene der Hydromechanik, wenn A= u = 0. gemacht wird. Zu den allge- 
meinen Gl. (IT) und (2) treten nun die Bedingungen hinzu, die die Lösungen als Rayleigh- 
wellen charakterisieren sollen: die Teilchenverrückungen nach den drei Achsen u, v, w sollen 
gegen das Innere des Mediums hin, also mit wachsendem 2, nach e-Potenzen abnehmen. 
Partikulare Lösungen dieser Art gibt es tatsächlich unter bestimmten Einschränkungen 
für A/, w. Wir bezeichnen sie auch hier, wie im rein elastischem Fall als freie Rayleigh- 
wellen [2]. 

Unsere Fragestellung soll nun so lauten: Lassen sich die einzelnen Partikular- 
lösungen derart superponieren, daß gewisse Anfangsbedingungen (und zwar 
solche, welche die Anfangslage, und solche, welche die Anfangsgeschwindig- 
keit der einzelnen Massenpunkte festlegen) erfüllt werden? Mit anderen Worten: 
Wenn zu Anfang der Zeitrechnung (t=0) in jedem Raumpunkt x, %, z die Teilchenverschie- 
bungen u, v, w sowohl wie die Geschwindigkeitskomponenten bekannt sind, kann es dazu 
Rayleighwellen geben, die diese Anfangsbedingungen erfüllen? Wenn dies nicht der 
Fall ist: wie viele und welche Anfangsbedingungen kann man vorschreiben, um gerade 
jenen Grad von Willkür auszuschöpfen, der in den Konstanten der Partikularlösungen 
enthalten ist ? 

Es ist ersichtlich, daß die ganze Frage von ähnlicher Art ist wie jene, die in der Theorie der 
schwingenden Saite behandelt wird. Die hier zu lösende Anfangswertaufgabe ist aber insofern 
schwieriger, als man es nicht mit einer, sondern mit drei Komponenten der Verschiebung’ 
zu tun hat, die nicht unabhängig voneinander sind. 

Bei der Lösung der Aufgabe geht man so vor, daß man, ausgehend von passend gewählten 
Partikularlösungen, diese durch Superposition zu allgemeinen Lösungen zusammenfügt und 
sie dann den Anfangsbedingungen anpaßt. Dabei zeigt es sich, daß der dreidimensionale Fall 
ein wenig anders zu behandeln ist als der zweidimensionale. 


1) Über die verschiedenen Möglichkeiten für die Wahl des Spannungstensors und ihre Berechti i i 

( gung vgl. die Arbeit [1] 
des Verfassers. Dort finden sich auch Hinweise auf die Literatur. Die Zahlen in eckigen Klamme i 
Schrifttumverzeichnis am Schluß der Arbeit. x x a ne ne 
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Das zweidimensionale Problem. 


1. Die Partikularlösungen. Die zu integrierende Differentialgleichung ist wieder {B), 
doch bedeutet 3 jetzt den zweidimensionalen Verschiebungsvektor mit den Komponenten 


u, w, bezüglich der x- und 2-Achse. Der Spannungstensor hat nunmehr bloß die drei Kom- 


ponenten n 
a RER 
= AO+2nEe + Ei O+2u a RER 
BRD, (3), 
o9w du ‚o (dw ou yet 
Prz a 1 +4 4 - 


so daß (I) zu integrieren ist unter der Annahme, daß | 
BEE ee ed BR: 


Aus den zahlreichen möglichen Lösungen von (I) greifen wir die Oberflächenplanwellen 
(Rayleighwellen) dadurch heraus, daß wir in dem.sonst üblichen Ansatz 


od 0W a 4 
a ee De 


die neu eingeführten Funktionen ® und Y mit z nach e-Potenzen abklingen lassen: 
DB — Ae- net it mid), W_ Be-n2+ i(a,©+ pri) 


Dementsprechend sind die Konstanten g, und g, reell und positiv vorauszusetzen. Was die 
anderen als willkürlich eingeführten Konstanten A, B, a,, a,, ?,, pP, anlangt, so zeigt man 
leicht mit Hilfe der Oberflächenbedingungen (4), daß a, =a,=a und p,=p,;=p ausfallen 
muß, wo nun aeine reelle und p eine im allgemeinen komplexe Zahl bedeuten, die wir ae 
in der Form 


p=p+tie p,erell).... rennen (7) 


voraussetzen. Ferner ergibt das Eintragen des Ansatzes (5), (6) in (I), daß A, u, 4, w nicht 
voneinander unabhängig sein können, sondern daß zwischen ihnen die Bindung besteht 


“=604, w=op (oa—=willkürliche Konstante) . ».. .%.. (8). 


In dem betrachteten Medium sind daher Rayleighwellen nur möglich, wenn 
der Deformationsvorgang überlagert wird von einem Fließvorganß, der durch 
die Größe o charakterisiert wird. Die Konstanten A’, w' charakterisieren demnach das Medium 
als ein „‚fastelastisches‘“, d.h. ein solches, in dem elastische Deformation und dazu propor- 
tionaler Fließvorgang nebeneinander gehen. Was aber A, B anlangt, so bestimmen sich diese 
beiden Größen bis auf einen willkürlichen Faktor genau mit Hilfe der beiden Oberflächen- 
bedingungen (4), die mit einem linear-homogenen Gleichungssystem für A und B gleichbedeutend 
sind. Da dessen Determinante verschwinden muß, erhält man eine zusätzliche Gleichung, die 
(unter der für die Erde weitgehend richtig gefundenen Annahme A —=u) auf die bekannte Gleichung 
dritten Grades 


3—1222+14x—4=0 ....:... BE (9) 
für jene Größe x hinausläuft, die in den Ausdrücken 
ar -"- +2 Lasse RE (10) 


- für die Dämpfung e und die Wellengeschwindigkeit V eine Rolle it Von den drei (reellen) 


Wurzeln dieser Gleichung ist nur die kleinste, «—=1—1//3 —0,4226 . .. brauchbar, denn nur 


sie macht die Exponenten 


2 ; ss 
en %=+tayl—2x EN EN ae) (11) 


reell. Unter Einführung der Abkürzungen \ 
eh? — (1 —x)er?—=T, (a2), — (1—x)e=h? Ler@?—1],(a2)... (12) 
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und einer willkürlichen Konstanten C kann man die beiden sich ergebenden Partikularlösungen 
unseres Problems in folgender Form anschreiben; Pr 


u=--CT; (a; +) e=*!sina(c+ Vi), 


nat (a, z2)e=*tcosa(x+-Vt), 
—2% 
oder 


u= +CT,(a,2)e”*tcosa(z +71), 


et 2) e-*tsin a(z + Vt) 
% 


Die gesuchten, freien Rayleighwellen erweisen sich mithin als gedämpfte Schwingungen, 
wie es bei der Beschaffenheit des Mediums als eines nicht vollkommen elastischen von vorn- 
herein zu erwarten war. Die physikalische Bedeutung der außer C noch auftretenden will- 
kürlichen Konstanten «a ist durch die Beziehung 


festgelegt, wo L die Wellenlänge bedeutet. Was die Spannungen betrifft, so gehört zur oberen 
Partikularlösung der Spannungstensor ; 


Prs = Dan ol ( +2) e=hı2 + (1— x) eanrleo« cosa(s+ Vi), 


! V a is 
2 = 2aul 1 re mern) erttsina(wt Vi) 1 
9%; =2anC (1 —x%) (eT7 FE N) e=*tcosa(s+ Vi) 
zür unteren hingegen der Tensor 
Pr = 2aul - (1-3): +(1— x) e7%2| e7*tsina(2 Vi), 
a), 


"leB2: = 2au co) 1- nl are) er ’tcosa(2+V!t), 


Pız =2aul (1 —x) (eT H?—e7%2) e=*tsina(2 +V it) 


Daraus ist unmittelbar zu ersehen, daß an der Oberfläche z=0 die Komponenten Pxz, Pzz 
verschwinden. 


2. Die Anf angswertaufgabe. Wegen der Linearität der Differentialgleichungen läßt 
sich die allgemeine Lösung U(x,2,t), W(x,z,t) aus den partikulären durch Superposition 
herstellen: / 

Merz) 2 (2er 

0 


[C,sina(& +Vt) +0,cosa(x +Vt)+0;sin a(® —Vt)+0,cosa(x—Vt)]da, 


S (17). 
W(x,z, 1) = KfT,(a, Ay gTmalg He \ 
) 
[C, cosa(& +FYt)— O,sina(x +Vt) +0,cos a («— Vi) — Q;,sin a(2c—V t)]da 
Hier wurde zur Abkürzung 
1—x au 

u a 1 ER 
\1—2x% | a 

gesetzt. Die Konstanten O,, O,,.... sind nunmehr als Funktionen von a anzusehen, Die Inte- 


grationsvariable a steckt ferner vermöge (10) und (18) in a und V, in V steht sie unter dem 
Wurzelzeichen. Wegen des Zusammenhanges (14) zwischen aund Z bedeutet dies eine Dispersion, 
deren Größe (außer von u und o) von der Konstanten o abhängt, Da diese aber sehr klein aus-' 
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fällt [1], so klein, daß für die in der Seismik auftretenden Wellenlängen eine Dispersion unter 
der Meßgrenze liegt, kann man für die Geschwindigkeit in erster Annäherung 


| B se 
VE RE Ka Aa a 


setzen. ‘Mathematisch bedeutet es eine wesentliche Erleichterung, die Integrationsvariable 
nicht unter dem Wurzelzeichen zu wissen. Eine weitere Vereinfachung erzielt man, wenn man 
fordert, daß sich alle Daten, die gegebenen sowohl wie die gefragten, auf die Oberfläche/2=0 
beziehen sollen. Die in (12) definierten Ausdrücke I‘; (a,2); I',(a,2) enthalten nämlich wegen (11) 
das a im Exponenten. Beschränkt man sich jedoch auf die Oberfläche, so reduziert sich jeder 
der beiden Ausdrücke auf die Konstante x, die man sich sofort mit den willkürlichen Funktionen 

 Cy(a) vereinigt denken kann. Unter diesen vereinfachenden Annahmen kann man die allge- 

:meine Lösung unseres Problems in der Form ansetzen 


5 | 6 U (s, t) f e-#at (cos ax [(C,— O,)sin aVt—(C, +0,) cos aVt) 
0} 4 
— sin ax [(0,— O,)sinaVt +(0; +C,)cosaVt]}da, 


a BR 
W(z,t)= Kfe- «e*t fcos aw |(O,— O,) sin aVt +(0,-0,) cos aVt] 
ö 


\ +sin ax [(O,— O,)sin aVt— (0, +0,) cos aVt]} da 
Die willkürlichen Funktionen C,, O3, ... bzw. deren Linearkombinationen 0, + (0, (+ C} 
sind nun den Anfangsbedingungen anzupassen. KEN ; 
Als Anfangswertaufgabe bietet sich zwanglos die folgende dar: N 
1. Im Zeitpunkt t=0 soll für die Vertikalverschiebung an der Oberfläche ein ‚,‚Profil‘ 
(x) vorgeschrieben sein - 
t=0, RE ae 0 lo Re (21). 


2. Die Verschiebungen sollen aus einer bestimmten Verteilung der vertikalen Anfangs- 
geschwindigkeiten g(x) hervorgehen 


0, )_,-@ A RE RE SEE NEN 


d.h. es ist für den Anfang.der Bewegung auch ein bestimmtes ‚‚Geschwindigkeitsprofil‘ für 
die Vertikalkomponente vorgeschrieben. 

Es zeigt sich, daß durch diese beiden Anfangsbedingungen die willkürlichen Funktionen C‘, (a) 
vollkommen festgelegt sind, insbesondere ist durch sie auch die Horizontalverschiebung U (z, t) 
zu jedem Zeitpunkt mitbestimmt. Die angeführte Anfangswertaufgabe ist übrigens nicht die 
einzig mögliche. Man könnte ebensogut die Horizontalkomponente statt der Vertikalkompo- 
nenten auszeichnen oder sich irgendeine andere, passende Zusammenstellung zurechtlegen. 

Aus (20) folgt für i=0 


fi) =[IK (0, +0,)cosac—K(C, +C,)sin ax] da, 
0 
& | ERBPEAE 
\ (2) = [S[K(C,— 0) aV-cosaw +K (0, —0C,)aV-sinax]da 
. 0 
—o<r<+ 


Um die Funktionen C‘,(a) so zu bestimmen, daß sie diesen Bedingungen genügen, greifen wir 
auf den Fourierschen Integralsatz in der Form 
© + 


oo +» 
"@-[lt] (McosahaA-oosaw +1. | Pia) sin aA dA-sinaz da... (24) 
4 
0 a — © 


zurück?). Identifiziert man hier F(x) das eine Mal mit f(x) aus (23), das andere Mal mit g(x), 


} 2) Die Voraussetzungen, die über die Funktion F(x) gemacht werden müssen, um die Existenz des Fourierschen 
Integralsatzes zu sichern, sind für die in der Geophysik auftretenden Funktionen in der Regel erfüllt. Um Sprungstellen 


der Form 
re +0) + F(&—0)] 


zu vermeiden, verlangen wir zusätzlich die Stetigkeit der Funktion F (x) im Intervall 0< x <». 


’ 
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so erkennt man, daß 


+ +» 
1 } 
+01 = | 1) os aadA, C; +0,= | fa sin «2.42, 
Br 2 Honst2ny 


+. +0 
1 i 1 
0,0 = gps Wsinarda, 0 G= zugp|oWeosarar 


zu setzen ist, Dadurch sind die Funktionen C,(a) so festgelegt, daß sie, in (20) einge- 
setzt, die allgemeine Lösung den Anfangsbedingungen anpassen. Man erhält als allgemeine 
Lösung ‚ 


x 


ge - . Be 
v0 |raaal [e-e=tsin aa V) da + |e-="sin UBER vi) da! 
9 


= .00 0 


to oo ® 
et (4) aa e-@e*tsina(A—x) sin aVtda, 
: (26). 


oo 


+ © Ex h 
W”ai)= 1 Jia) dA Lfe-wetcos a(A— x—Vt)da +[e=e@toos a(i— x+Vt) dal 


—o 2 


() 
Were 
ER, I p— qa:att er 3 da 
N +, |swar|te cosa(A— x)sinaVt da 
; —® 0 


Die hier auftretenden Integrale sind nur zum Teil elementar auswertbar. Man findet, wie nicht 
weiter ausgeführt werden soll, y 


oo 


1 _a—aH+ vo: 
geh Ts F} 
J.= — adatt Niet N 4 a®t 
= IE cosa(A— xz+TVt)da >77 Mar ; 
0 
1— ak vi Ser Yale 
2aVt 


_— 
x 


1 Er 
el = [we 1 ARTE Een De‘ 4a:rt +yi 
9 e sina(A— x +TVt)da arı® e+t’’dy 
0 0 


Während J, (das übrigens in mancherlei theoretischen Untersuchungen auftritt) in dieser Form 


belassen werden kann, ist es bei J, zweckmäßig, sich einer Reihenentwicklung zu bedienen. 
Ferner läßt sich zeigen, daß 


ı— a— Vi 
es "2avi 
= | "et cosaa—asinaVvtda= Et | e-ray. BEE v7) 
x | 


A—a+Vt 
2aVt 


ausfallen muß, während 


oo 


SR AR 
= [Le eetsin a(A— x)sinaVtda 
0 


sich nicht weiter auf eine einfachere Formel zurückführen läßt. Auch hier wird man zu einer 
Reihenentwicklung greifen müssen. Was J, anlangt, so ist es stets als Summe von Gauß- 
schen -Fehlerintegralen darstellbar, die bekanntlich tabelliert vorliegen (vgl. z.B. Hayashi: 


Fünfstellige Funktionentafeln, Berlin 1930). Die Lösung der gestellten Anfangswertaufgabe 
sieht daher so aus: 


. 


a 


> BE ih 


ae 


BEE Ss 
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2 
z oo 


Sb @ tm ‘ ar 
Va=5.z | SWmarle-wetsinaa— a VHda + fe #stsinaa— = F Pi) da} 
0 » h 


io 0 N 


— DD 


+% = ' 
Hop | mar [Le #etsin aaa) sin aVtda, 
0 


A—ı—Vt 


I  adszvg Maar EL di 

ee aan He u +4 | owar |eray 
daynt 272 Es 
Gr an A— ı+Vt 


| t NW : 
Die in dieser Lösung auftretenden Funktionen, die zunächst ganz allgemein sein können, müssen 
nun, nach Maßgabe-der gerade vorliegenden Aufgabe, entsprechend spezialisiert werden, Dabei ° 
ist aber zu beachten, daß sie so gewählt werden müssen, daß die entsprechenden Integrale 
tatsächlich existieren. Eine Möglichkeit der Spezialisierung. besteht z.B. darin, für i=0 zwar 
ein bestimmtes Profil für die Vertikalverschiebung an der Oberfläche vorzugeben, die Bewegung 
des Mediums aber aus der Ruhe heraus erfolgen zu lassen. Dann hätte man g(x) = zu setzen _ 
und (30) würde sich wesentlich vereinfachen. Eine andere Möglichkeit bestünde darin, anzu- 
nehmen, daß für i=0 keine Vertikalverschiebung vorliegt, also f(x) =0 wird, daß aber 
die Oberflächenteilchen in diesem Augenblick einen vertikalen Impuls bekommen, der ihnen 
ein Geschwindigkeitsprofil g (x) vorschreibt. Auch in diesem Falle würden sich die Formeln (30) 
vereinfachen, besonders aber dann, wenn man von dem Impuls voraussetzt, daß er nur in einer 
kleinen Umgebung des Ursprungs erfolgt, die man dann noch beliebig verengen kann. Dieser Fall, 
die. punktförmige und symmetrische Anregung des Mediums, bietet sich in der Praxis sicher dar, 
allerdings erst bei der Behandlung des dreidimensionalen Problems. Das Fallen eines schweren 
Körpers auf. die Erdoberfläche, eine Explosion oder ein Bombeneinschlag liefern Beispiele dafür. 

Eine weitere Vereinfachung in den Formeln von (30) läßt sich erreichen, wenn man zum 
Grenzfall des vollkommen elastischen Körpers übergeht. Dann hat man lim e—0, zufolge (18) 
mithin lim a —0 zu setzen. Die Lösung (30) wird.dann zu 


+ © © co oo 
Ua) = EN (A) aa [sin aaa Vidatz.g [fa [sin aller Vida 
= {oJ } _—» 0) 


1er fi 
urn [war] sin a(A—x)sinaVtda, 
1 +0 oo 1 + Be (31) 
Win, |swar|enaa— 3 Püdatz,|smaa|saa—3F Vt)da 
— © 0} —o» 
a A 
+apfowmaa|[ cos a(A— x)sinaVt)da 

7 ga 
Liest man den Fourierschen Integralsatz (4) von rechts nach links, so. erkennt man, daß 
die erste Zeile im Ausdruck für W zu 

- 1 - 

5 f@—V)+f(@+ PN) 


wird, während die zweite Zeile noch einer geringfügigen Umformung bedarf. Es ist nämlich 


0 füry>l1l, 

7 
A rn u ” y=1, 

1—4 

[7 sin avı-e0s aa— 2)da= | sin acos ayda= AR 1>y>—1,wy=-—— (32). 
Ja ya 2 Vt 
ö 

7 

Da v=r-l, 

0 ” Vz 
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eines jener Integrale, die von einem Parameter (hier y) unstetig abhängen?). Im Hinblick 

auf die Abhängigkeit des Parameters von x und ? kann man auch sagen, das Integral habe 

den Wert E BR 

3 ak | ER 0 
2 4 Le, 4 | 

A>z+Vi, i=2+Vi, = —-Vi<i<e+Vt, A=s2—Vi, A<s—Vi. 


Zum Doppelintegral der zweiten Zeile liefern daher nur jene Werte des A einen Beitrag,_die 
zwischen den Grenzen &— Vt und x + Vt liegen. Es ist somit f 


+» 5 Bi vt 
[owar| z eosaa—a)sinavı da= D2 g(A)dA., 
aa 0 s—Vt 


. Das entsprechende Doppelintegral im Ausdruck für U läßt sich ebenfalls auf ein einfaches 
Integral zurückführen. Es läßt sich zeigen, daß 


Be RR Are TER 
fasn« tsin a(A— x) a=7log Enz 
ausfallen muß. Für die Verschiebungen U (z,t), W (z,t) des Bodenteilchens findet man daher 


+» 20 © 
rad I [e-@ersin a(A— x — Vi) da+ fe «tsin a(A— 24V) dal 
Re 


ER [) 


+® N. 

1 Ve ee: 
tage % a az, 

imRv |" ie = 


(33). 


. 


1 s+Vt 
Wange) ++ VO) + gr [oa 


zs—Vt 


Dies ist die Lösung des Anfangswertproblems für den Fall, daß das Medium 
als vollkommen elastisch vorausgesetzt wird. Der Ausdruck für W ist deshalb be- 
merkenswert, weil es derselbe ist, der auch in der Theorie der schwingenden Saite auftritt, 
wenn man das entsprechende Anfangswertproblem löst (Saite unbegrenzt, gegeben Anfangs- 
lage und Anfangsgeschwindigkeit der einzelnen Saitenpunkte). Die Deutung der Lösung ist 
hier wie dort gleich, so daß es sich erübrigt, darauf einzugehen. Erfolgt die Bewegung aus 
der Ruhelage heraus, so ist f(x) =0 zu setzen und die Lösung wird zu 


ee en s+ Vt 
1 A— avi: N 
Ua) = j 10 1 | BDRNEN-t 
were |, an PET (34) 


Das dreidimensionale Problem. 
3. Die Partikularlösungen. Zu integrieren ist jetzt (I) mit den Oberflächenbedingungen (2). 
Der Ansatz 


oD* er 


0x 92092)’ 


u=eirl 


ven?" ne ung. [aB" X 


EERTIN en we: 2 

oy Hayes 02 Ar rk x) le, 
führt dazu, von den neu eingeführten Funktionen ®* und X, die nunmehr von der Zeit nicht 
mehr abhängen, das Erfülltsein von 


Ad» FTMROr=D, . AXEEK ED Rn a . (36) 
zu verlangen, wo 
a is ER 4 ale a 
Ian u) k KT, DD TIEREN 


°®) Wegen der Behandlung dieses Integrals siehe z.B.B. Riemann: Partielle Differenti i 
£ «PB. BD: : ntialgleich 5 - 
ausgabe von Fritz Emde, Braunschweig 1938). a 
% 


TEE WERE 
PAR 


° Koordinatenachsen, sondern die Radiaiverschiebung u; und die 
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> * x 


gesetzt wurde. Da die Nerschiöhunden u, v, w zu Oberflächenwellen gehören sollen, wird man 
D+= e-t2K(v}y), X=et%:K (x, y) ’ 


setzen. Die Konstanten g, und g, sind wieder reell und positiv vorauszusetzen.. Zusammen 


mit den Se NS erhält man En für X (x, y) die folgende Differentialgleichung 


..»92K 
ae} ar Fate E= On un Pla deaer Reel 
e®=g Yon, Me RENT Se u N: 


Der ze zu Polar- bzw. zu Zylinderkoardinaten vermögee=rcosd, y=rsind,2=2, 
ji ‚tzi + ® bewirkt, daß sie in 
PR LOK RK 
a? 
Er Er Or Aurrr) a ir 


uherzent: Wie beim Problem der schwingenden Membran setzt man nun 


cosn® 
K=Jn(ar) BERN n>(0, ganz, 
‚und erreicht Burch die Forderung, daß n ganzzahlig sei, daß K eine eindeutige Funktion von » 
wird: bei Vermehrung des 9 um 2 reproduziert sich der Wert der, Funktion‘). 
Für J„(ar) aber ergibt sich als Forderung das Bestehen der Besselschen Differential- 
gleichung 
d? Jn 1 dm | n?2 a 
— + —— - — )M"=0... 2. 40), 
d(ar)? EHRT (ar) et (ar)? ne h Ser 


wodurch sie als Besselsche Funktion vom Index » definiert ist, und zwar als eine solche erster 
Art, denn die Besselschen Funktionen 'zweiter Art oder Neumannschen Zylinderfunktionen 
genügen zwar auch dieser Differentialgleichung, sind aber physikalisch unbrauchbar, da sie 
an der Stelle »—=0 logarithmisch unendlich werden. Unter Ans Bn, An, Bn willkürliche Kon- 
stanten verstanden, wird mithin 
87, —e7-ı2(Ancosnd + Bnsin nd) Jn(ar), . 
=. 021,2 7,: en a eh RRBAAFE TEE 
X=e71%2(A,cosnd+ B,sinnd) Jn (ar) 
und. damit sind die Verschiebungskomponenten , v, w fest- 
gelegt. Nach Einführung der Zylinderkoordinaten interessieren 
aber nicht mehr die Verschiebungskomponenten längs der 


Tangentialverschiebung us, die, wie man sich an Hand von 
Bild 1 leicht überzeugt, mit u, v, w durch 


Y,= Ucosd+vsind, 
UY=—-usin®-+vcos” 


zusammenhängen. Ebenso wird man nun nicht mehr nach den 
Spannungskomponenten Pzz, Pzy, . . . fragen, sondern nach den 
Komponenten Prr, Prs, P68, Prz» Paz Prz. Die elastische Span- 
nung hat den Charakter einer Tensordichte, d.h. die Transfor- 
mationsformeln für die einzelnen Komponenten beim Übergang 
von einem (beliebigen, auch krummlinigen) Koordinatensystem | 
zum andern sind mit der Funktionaldeterminante der Transfor- 
mation behaftet (auf diesen Umstand hat erstmalig L. Bril- 
louin [3] hingewiesen). Bei Berücksichtigung dieses Umstandes 
ist es nicht schwer, die fraglichen Spannungskomponenten auch tatsächlich zu berechnen. Im 
Gesamten lautet das Ergebnis: 


Bild 1. 


% Verschiebungen 
U, = aeipt (Mne- 42 — g, Mye-%2) J,(ar), | 


N = — eint (Nne-n2—q, N;e— 42) In (ar), . (43). 
w= eirt[— q, Mne=-u2-+ (9, + R?) Mne=%2] In (ar) | e 


*) An sich wären auch negative, ganzzahlige Werte von n zuzulassen, doch genügt im folgenden die Berücksich- 
tigung der n 20. 


RT EN REES WENN AUF BRNGE NER ENE 


N. 
N an 


= . 


Fr 


* 


Spannungen | 
Dr: = wart [— 2g, Meran + 2, +8) M',e-%2] J,, 


ps=pten ng New t+@gHRNield |...) 


Hier wurden die Abkürzungen BR a : 
Ancosn® + Bnsinn® = Mn, — Ansinnd + Bncosn® = N, ee 
Ay,cosn®-+ B,sinn®= M}, — A,sinnd+ Bncosnd—=Nn 


verwendet. Die weiteren Überlegungen verlaufen wie im zweidimensionalen Fall. Aus den 
Oberflächenbedingungen läßt sich, bis auf einen willkürlichen Faktor genau, das System der 
Konstanten An, Bn, A}, Bj} berechnen. Das Verschwinden der Determinante des entsprechenden 
Linearsystems führt auf dieselbe Gl. (9) für x, so daß man für die Dämpfung e und die Wellen- 
geschwindigkeit 7 wieder die Ausdrücke (10) erhält. Auch die Ausdrücke (11) für 41 und Ge 
bleiben unverändert. Die gesuchte, partikuläre Lösung des Problems nimmt bei Berücksichtigung 
all dieser Umstände sowie Zusammenziehungen die Gestalt an 


= ir alla,2)e=et[cosaVt( An cosnd + a,sin nd) | f 
—sinaVt( bn’/cosn® + bi sin nd)) J„(ar), 
U= - T', (a,2)e—*t [cos aVi(— sinn® + a) cos nd) 


—sin aVt (— bnsinnd + 5, cos nd)] Jn(ar), 
# T,(a,z)e —s[cosaVti( aAncosnd +a,sinnd) 


; —sinaVt( bncosnd +b,sinn®)] Jn(ar) 


Entsprechendes hätte man für die Spannungen durchzuführen. In der Tat gehören Verschie- 
bungen und Spannungen insofern zusammen, als zu jeder, durch Superposition gefundenen 
neuen Lösung auch neue Spannungen gehören, die man durch Superposition der. alten Parti- 
kularspannungen findet 5). 

Daß schon einer einzelnen Partukilarlösung aus (46) physikalische Bedeutung zukommen 
kann, ersieht man aus folgendem Sonderfall: denkt man sich”die Anfangsbedingungen so be- 
schaffen, daß das Erschütterungsfeld um den Herd zylindersymmetrisch ist (wie es etwa bei 


einer Sprengung an der Erdoberfläche, einem Bombenwurf, einem schweren Fall eintreten ° 


wird), so muß die gesuchte Lösung frei sein vom Azimut d. Wir erkennen aus (46), daß schon 
die nullte Partikularlösung zur Beschreibung solcher Vorgäng& herangezogen werden kann. 
Man erhält nämlich 


u= — I;(a,2)e-et(a,cosaVt—b,sinaVt)J,(ar),, w=0, 
w= SH T',(a,2)e—*t(a,cosaVt—b,sinaVt) J,(ar) FEN 
Y1—2x n 
Hier wurde auf die aus der Theorie der Zylinderfunktionen bekannte Beziehung [5] 
J,(ar) = —Jı(ar) 


zurückgegriffen und außerdem der sowohl in u, als auch in w auftretende Faktor x mit den 
willkürlichen Konstanten vereinigt“). Die in den Konstanten a,, b,, a liegende Willkür kann 


5) Was die Partikularlösungen (46) anlangt, so unterscheiden sie sich von den vom Verf. aufanderem Wege [4] ge- 
fundenen Lösungen 


u=fl'(a,2)e ®tcosa(y + Vt),  u=—flı(a.)e”sina(y + Pi), 
v=gl,(a,2)e*tcosa(y + Wit), v=—gT', (a,2)e *tsina(y + Pt), 
a(1—x) e 2 a(1— x) __ 
w= _- = T,(a,2)e ®t sina(y + Pt), w= IT, (a,e)e3°t 608 rt), 
02 ei ( ) er 2.(a,2)e”, a(y + Pt) 
3 y- n cz + 2. Y; a=Y/f:+g, f, go, willkürliche Konstanten, 


und es wäre der Nachweis zu erbringen, daß sie diesen nicht widersprechen. Dieser Nachweis gelingt, wenn man die letzt- 
genannten Lösungen, die Planwellen darstellen, durch Superposition zu Zylinderwellen zusammenfügt. Die im Text ge- 
gebenen Lösungen (46) erweisen sich dann als Folgen der hier angeführten. Vgl. dazu die Darstellung der Theorie der 
Zylinderfunktionen von R. Weyrich [5], in der der Versuch gemacht wird, die Erzeugung von Zylinderwellen durch 
Überlagerung von Planwellen zur Grundlage der Darstellung zu machen. ‚ 

*) Das Ergebnis ist in Übereinstimmung mit einem schon früher [4] vom Verfasser hergeleiteten, erscheint hier aber 
als Sonderfall der allgemeineren Partikularlösungen. 
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dazu Henke werden, eine zu suchende Lösung durch eu den aylindersymmetrischen 
Anfangsbedingungen anzupassen. 

4. Die Anfangswertaufgabe. Um die allgemeine Lösung U,8,2,1), Volr,d, 2), 
W(r,®,2,t) aufzustellen und dann beliebigen Anfangsbedingungen anzupassen, wird man die 
ganze Willkür ausnutzen, die in an, @,, Dr, b,, a enthalten ist, d.h. man wird alle Partikular- 
lösungen (47) überlagern. Dies bedeutet, daß man nicht nur die zu den verschiedenen Indizes 
n=0(,1,2,... gehörigen Lösungen durch Reihenbildung superponiert, sondern daß man auch 
eine Überlagerung der verschiedenen, zu a gehörigen Lösungen, also eine Integration nach a 
vornimmt. Darnach ist es leicht, die allgemeine Lösung hinzuschreiben. Wie im zweidimen- _ 
sionalen Fall machen wir wieder vereinfachende Annahmen: Gegeben sowohl wie gefragt sollen 
nur Daten von der Oberfläche z=0 sein, die Wellenausbreitung sei dispersionsfrei [es gelte 
also (19) statt (10)]. Die so verstandene Lösung werde mit U,(r,9,t), Us(r,d,t), W(r,®,t) 
bezeichnet und sieht Ar ke aus: . 


U,(r,d,t) = far« 3 [cos aVt( er + a,sin n®) 
"= 0 _ sin aVt( bncosnd +b5,sinn®)] J,(ar)da, 


Vo(r,9,1) = >, [+10 aVt(—a„sinnd® +a,„cosnd) er ARTE 
a — sin aVt (— bnsinn® +, cos nd)] Jn(ar)da, 


W(r,®,i) = xf ae—et 53 [cosaVt( a„cosn® +a,sinn® 
"= sinaVi( bncosn® +b,sin N) 
Wir lösen folgende Anfangswertaufgabe: 
1. Im Zeitpunkt i=0 soll für die 'Vertikalverschiebung an der Oberfläche ein ‚Profil‘ 
f(r,®) vorgeschrieben sein: 


iQ, TERN EN MER MARIN SEN (49). 
2. Die Verschiebungen sollen aus einer bestimmten Verteilung der vertikalen Anfangs- 
geschwindigkeiten g(r,®) hervorgehen- 


&, Ban REIT NE NEE (50), 
ö1 t=0 
d.h. es ist für den Anfang der Bewegung auch ein bestimmtes, an die "Oberfläche 
gebundenes ‚‚Geschwindigkeitsprofil‘‘ der Vertikalkomponente vorgeschrieben. 
Gefragt ist nach der Teilchenverschiebung zu einem beliebigen Zeitpunkt. Auch hier, beim 
dreidimensionalen Problem zeigtessich, daß durch diese beiden Anfangsbedingungen die willkür- 
lichen Funktionen vollkommen festgelegt werden können, daß also die übrigen Komponenten nicht 
mehr freiverfügbarsind. Ist der vertikale Verzerrungszustand festgelegt, so sind dadurch auch die 
horizontalen Verzerrungen mitbestimmt und mit ihnen auch die Spannungen, die ihren Ausgleich 
in den Schwingungen des Mediums finden, deren räumlichzeitlichen Ablauf eben die Lösung des 
Anfangswertproblems beschreibt. Dabei ist unter ,‚Schwingung‘ nicht nur eine oszillatorische Be- 
wegung zu verstehen —denn bei unbegrenztem Medium wird man von vornherein kein Oszillieren 
erwarten können —, sondern ganz allgemein die Teilchenbewegung unter dem Einfluß einer An- 
fangsanregung. Auf PAS angewendet, hat man es also mit folgenden Anfangsbedingungen zu tun: 


6,0 fa ( (An cos nd + sin nd) In(ar)d a, 


n=0 
(SE) 
9(r,9) = Se a > l(ean+aV bn) cosnd®-+(ea,+aVb,)sinnd] Jn(ar)da 
N) n=0 r 
Um nun die Funktionen a„(a), a,(a), ... so zu bestimmen, daß diese Bedingungen erfüllt 


werden, greifen wir auf den bekannten Satz aus der Theorie der Besselschen Funktionen 
zurück, der einerseits eine Verallgemeinerung der Fourierschen Reihe, andererseits eine solche 


. des Fourierschen Beh darstellt und Ra Ra geschrieben werden kann: 


1 
a F(r,ß)d 
F (r,®) Dan (ar) Jada|y. (ar) jrar | (t,P)dß 
© © © 2n 
1 5 
2 I ar) Jada | In (ar)rdr [Fo u, cosna 2... (02). 
u.) co 2n 
3 | Intarsada [ Intarrar | F6.M sinnßdß-sinn® 
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Die erste Gleichung von (51), die man auch in der Form anschreiben kann 
fd) = K[J,(ar)ada'a,(a) 
N) 


+ KS [J,(ar)adaa,(o)-cosnd+K [Intar)adaa, (a) sin nd 2 
n=10 ; n=1 
zeigt nun beim Vergleich mit (52) und Identifizierung von F(r,0) mit f(r,d), daB “ 
1 \ 5 1 e | 
An (a) = mr a, (a) = SEMu (a), n=0,1,2,... 


Ausfallen muß, wo nun 


1 = a Ta Pie | St. B)aß, 


Li (o) =/s In(ar)rdr ji firsB) 208 n BON N ne (53). 
a RR 


male) A fr.B) sinnßdß 


gesetzt wurde. Ebenso bestimmt man durch Entwicklung von g(r, 9) die Koeffizienten 
bn (a), du (a): 
1 1 


1 ra) (a) =— — gyMld), n=0,1,2... 


mit 
© 27 
Ei) = 2 [Istaairdr {te + ent map, 


0 


Lo) zn 
L2 (o) = [Untar)rdz FePß)tegs,ß)]eosnßaß,!i....... (54). 
== 12,30% 


Ma 2 fm anrdr] [f,ß) + eg(r, ß)]sin nßdß 


Damit ist die gestellte a gelöst, denn man hat nur die gefundenen Aus- 
drücke in (48) einzutragen. Das Ergebnis läßt sich in folgendem Satz ausdrücken: 

Es möge die Verschiebung eines, der Oberfläche z=0 des Mediums an- 
gehörenden Massenpunktes in Richtung des Radiusvektors mit U,(r,d,t), in 
Richtung senkrecht dazu mit Us(r,d,t) und in vertikaler Richtung mit W (r,d,t) 
bezeichnet werden. Wird nun die Vertikalkomponente im Zeitpunkt t=0 
den Bedingungen (49) und (50) unterworfen, so sind die Verschiebungen zu 
jedem anderen Zeitpunkt gegeben durch‘ 


U,„(r,d,t) el ae® J;(ar) LR (a) cosaVtda + f e=tJ4lar) DR (a) sin aWidalcos nd | 
J | 


n=0 0 


ns 


1 =» oo 1 x 
en = f ae: J;(ar) M) (a) cosaVtda-+ „| et], (ar) M) (a) sin aVtdaysin nd, 
n=0 0 


1 oo 
U;(r,d, = van ale stJ„(ar) M'% (a nal: alten (ar) MR X a)sin a Vidalcosnd 
co (55). 


1 
alle e! J„(ar) “(a )cosaVtda-+ une e®tJ„(ar) DE (a) sin aVtda!sinn®, 
0) 


SI e 


N 


1 
W(r,ö,t) = BI ae st J„(ar) 6% (a) cosaVtda-+ — f e-*t,J(ar) L% (a) sin aVtdal cosnd 


1 oo 
+2 1er st.J„(ar) M'% (a eos arida+ | et J„(ar) L% x (a) sin aVt da! sin nd 
n=0 0 


— en — —  — — — — — 


2. angew. Math. Mech. - | 
Ba. 25127 Nr.1 April‘ 1497  Bellich, Die Randbedingungen der Airyschen Spannungsfunktion 13 


Die hier auftretenden Funktionen L®(a), ... sind durch (53), (54) bestimmt. Bei zylinder- 
symmetrischen Anfangsbedingungen werden f(r,®) und g(r,®) frei von 9. Die Vereinfachungen, 
die sich dadurch ergeben, sind leicht ersichtlich. Die Lösung der Anfangswertaufgabe sieht 
so aus: er Ä hi; 


% 


U,(r,)= far, (ar) L®(a) cosaVtda+ ve (ar) L@ (a)sin a Pida! 
RT! : ö 


sK 


Uslr,t) =0, 


v 


(56). 


% 


D 


Wert)= a fee” "Isar 20 (a) cos aVtda-+ [een L@) (a) sin a ida! 
6 ö 


\ 


So wie die Anfangsbedingungen, ist auch die Lösung zylindersymmetrisch. Die Bewegung aller 
Massenpunkte gleicher Herddistanzerfolgtgleichartigund keiner verläßt, wie aus dem Verschwin- 
den der Komponente U;(r,t) hervorgeht, die Vertikalebene durch den Radiusvektor. Die Lösung 
(96) hätte übrigens auch unmittelbar aus (47) durch Superposition hergeleitet werden können. 
Was die Anwendung der hier abgeleiteten Resultate auf die Praxis anlangt, so könnte 
man etwa an folgende Fragestellung denken: Eine Vertikalverschiebung der Massenteilchen 
der Oberfläche zu Anfang der. Zeitrechnung ‚findet überhaupt nicht statt, d.h. es ist 
/(r,9) = f(r) =0. In einer bestimmten, kreisförmigen Umgebung des Ursprungs (des Herdes), 
die beliebig klein angenommen werden kann, erfahren die Massenpunkte in der Vertikalen 
einen plötzlichen Stoß (z. B. Fall einer schweren Kugel). Die Funktion g(r,9) wird dann nicht 
nur von r allein abhängen, sondern auch so beschaffen sein, daß sie nur innerhalb des genannten’ 
kreisförmigen ‚Bereiches nicht identisch verschwindet. In der entsprechenden Lösung 


U) = ge *!J,(ar) sin a Vida | Jotan)ota)rar, UHlrt)=0, 
/ l 
WAr,t)= i [er sen sin aVtda [Iotanlaterdr 
ö ö 


werden dann die oberen Grenzen der inneren Integrale auf einen endlichen Wert herabgesetzt 
werden können. Unter bestimmten Voraussetzungen über das Unendlichwerden der Funktion g(r) 
an der Stelle r=0 kann man weiter zu dem Grenzfall übergehen, daß der Radius des Kreis- 
bereiches unendlich klein, die angeregte Stelle also punktförmig wird. Die entsprechenden 
Lösungen der Anfangswertaufgabe werden dadurch weiter vereinfacht. 
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Die Randbedingungen der Airyschen Spannungsfunktion 


bei vorgegebenen Randverschiebungen. 
Von F. Rellich in Göttingen. 


1. Das Verhalten einer Platte, deren Durchbiegung.nicht mehr klein ist gegenüber der 
Plattendicke, wurde in den letzten Jahren mehrfach mit Hilfe der beidenvon Kärmän- 
schen Differentialgleichungen!) für die Spannungsfunktion ®(x, y) und die Durchbiegung 
w(x,y) beschrieben: 

AAB=E (wiy — Was Wyy) : 
ER? RE a 
12 (1— 2) AAw= Dyy war— 2 DayWwayt+ Dax wyy 


ı) Th.v. Kärmän, Enzyklopädie der mathematischen Wissenschaften Bd. IV Art. 27 insbesondere 8. 350 (1910). 
Literaturübersicht bis 1939 bei K. O. Friedrichs und J. J. Stoker, The nonlinear boundary value, problem, of 
the buckled plate, Nat. Acad. of Sciences, Vol. 25, S. 535—540 (1939). 
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(E Elastizitätsmodul, % Plattendicke, v Poissonsche Zahl). Die (über die Plattendicke ge- $2 4 
_ mittelten) Verschiebungen u (z, y) und v (x, y) ergeben sich, sobald ® und w nach (1) bestimmt se 
sind, aus den Gleichungen h | | 


MR a e 
y- Hy arm) Ne RER EAN 
2(14» | 
uUyt+ Ve — E N 


und zwar ergeben sich alle Lösungspaare von (2) aus einem bestimmten Paar u, v in der Form 
WLCgY+ 6» d!— C9% + C, mit geeigneten Konstanten c,, €, €;. Die (über die Plattendicke 2 
gemittelten) Spannungen sind 0, = Byy, oy= Da, T=— Dy 

Die Gleichungen der linearisierten Theorie entstehen, wenn in (1) und (2) die in w qua- 
dratischen Glieder gestrichen werden. 

Als Randbedingungen treten zu (1) zunächst die Bedinungen für die Durchbiegung w 
je nach den am Rande der Platte vorgeschriebenen Verhältnissen (eingespannte Platte, frei 
gelagerte Platte, usw.). Außerdem sind am Rande der Platte entweder die Oberflächenspan- 
nungen 9x, ?y oder die Verschiebungen a,» vorgeschrieben. 

Im Falle gegebener Oberflächenspannungen bedeutet das bekanntlich, daß ® und seine 
Normalableitung ®„ auf dem Rande vorgeschrieben sind, und zwar wird, wenn der Rand 
aus einer geschlossenen Kurve besteht, im Randpunkt, der zur Bogenlänge s gehört?) 


= / My) —yY(8’)) Pe ()— (el) — x (8) Pr ()]d Harz +by+e 


4 8 8 
D, =] Pr ds! — &n | By ds’ +axm + byn 
) 
Die Konstanten a, b, c bleiben willkürlich, können aber durch die weiteren Bedingungen 
If Bazay=0 | J[ 9242 dy= 0 If ®rdady=0 


eindeutig festgelegt werden; K ist das zugrunde gelegte &, y-Gebiet (Platte im unbelasteten 
Zustand). Es ergibt sich also hinsichtlich der Randbedingungen für Ö und w kein Unter- 
schied gegenüber der linearisierten Theorie. \ 

2. Im Falle vorgeschriebener Randverschiebungen?) %, v ist das anders, und es sollen 
im folgenden die Randbedingungen angegeben werden, die für Ö und w zu stellen sind, damit 
die nach (2) errechneten Verschiebungen u(%,y), v(z,y) am Rande die vorgeschriebenen 
Werte annehmen. Diese Randbedingungen sind, wie mir scheint, auch für die linearen Glei- 
chungen neu und geben dann an, welche Bedingungen beim ebenen Spannungszustand von 
der Airyschen Spannungsfunktion am Rande erfüllt sein müssen, wenn die Verschiebungen 
am Rande vorgeschrieben sind. Es seien also ©, w, vw, v vier Funktionen, die den Gl. (2) ge- 
nügen. Der Rand des zugrunde gelegten Gebietes K in der x, y-Ebene sei gegeben durch 
eine geschlossene Kurve = x (8), y=y(s), wo s die Bogenlänge des Randes bedeutet; 
Con = Ys» Yn = — %s Seien die Komponenten der nach außen weisenden Randnormalen. 
r«' Der Gedankengang der folgenden Rechnung ist dieser. Mit vw und v sind auf dem Rande 
auch die Ableitungen u,, v; nach der Bogenlänge gegeben. Kann man also u,, v; auf Grund , 
der Gl. (2) allein durch ® und w ausdrücken, so hat man für ® und w notwendige Rand- 
bedingungen gefunden. Der Symmetrie halber wird das nicht für ,, v; sondern für die beiden 
Kombinationen us &s + vs ys und Us &n + ds Yn durchgeführt. Die erste läßt sich unmittel- 
bar durch ® und w allein ausdrücken, die zweite erst nach Differentiation nach s. 

Es ist 

Ust + YY—= Us; + (Uy-+ ve) % y + vy Y; 

also nach (2) auch 


H v 1 
Us %s + vs Yy— M (Dyy % —2 Days Y-+ Pary) — (Deu + 2 Day %s Ya+ Dyy Y)— 5%; a: 


a N 0 S y 1 
I! (Dra &n + 2D2y n Yn + Dyy Yyı) ei) (Dirt Pu) — ZW; A 


?) Vgl. z.B. Biezeno-Grammel: Technische Dynamik, S. 99—100 Berlin, 1939. 

®) Vorgeschriebene Randverschiebungen in den einschlägigen Untersuchungen von K. M arguerre, vgl. insbe- 
sondere K.Marguerre und A. Kromm: Verhalten eines von Schub- und Druckkräften beanspruchten Platten- 
streifens oberhalb der Beulgrenze. Luftf.-Forschg. Bd. 14, S. 627—639 (1937) und die dort angeführte Literatur. 
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Bezeichnet # die Krümmung der Randkurve, also 


% = — Hm = + Hy Yıs = — Hyn = Ms, 
so wird - / 
Dis = — (Dorn + 2 Dayn yn + Du yn) + AB —#Dn, 
also 
1 ; 1 
ut +0 40-77 BE (3) 


Das ist eine erste Beziehung, der ® und w genügen müssen, wenn %, v auf dem Rande vor- 
geschrieben sind. 


Weiter wird 
Us&n + Vs Yn = Ur %s Un + UyYs in + Vz %sYn + %y Ys Yn 


’ % 
= — Ur Ys Yn -H (w-+ At Aue 


1 
lt nn (Pas E%s + Day Yan % + 2 yn) + a YnYı) . 


1 
— dy%g%n + 2 (Uy — dr) 


1 (ut Ys + wawy (Y 2) a y)+5 (m yv— da). 
Es ist 


(Dy)s = (®: In C D, Yn)s 5 Dr Un %s # o, (X, Y%; Fr %; Yn) By Day Yn Ys + ®; Yss —D, Vgs 
und daher 


5 Drain% + Dyy (in Ys + %s Yn) + Dyy Yn Ys = (On)s —x9;. f 
1+v 
E 


Also 


1 x 2 1 
Un + Y mn = — (Bu), — 9.) — Do (wz %Ys + Wz wy (Y — %s) — Wy 2%) +3 (Uy— da). 


Liegt.für s = s, eine Ecke vor, so ist jedenfalls 


s=s,+0 1 Er 1 s=8,+0 
in tt; pm)" = er (Da); —xD,). ar (wnw) EN )) 
re in a an 


eine zweite Bedingung für ©, w auf dem Rande. Eine dritte Be sich durch Differentiation 
der vorletzten Formelzeile nach der Bogenlänge: 


1+» £ 
(u zn + u = (Or) — (Di) 
= 
eE >, (Ugy%s + UyyYs— Yon %s — Yzy %) 
1 2 2 B 1 2 SUB 2 
ga 2 [2w, U; Ys Sr w, (X; Ya %;)] a} Wy [2w, X Ys Sr W, (Y— %s Ys) + 
+, (, 9 — %)—2w,% Y,] 


BR] Yyy [w, ( — 1 Yı) 2 W, %; Ys] 


1 
m [w; (Hs Y En %s Yss) iR 2w, W, (Y Ys —% Vs) — uy (Lips Y; ni %s Yss)] 


Auf der rechten Seite lassen sich die Ableitungen von u, v wegen (2) herausschaffen: 


(U, &n + 9%, Yn)s = ee. (Da) —rD,)s+ an (Dyyy le — 2 Payu Y— Pazz Y%t 2Dyay Ist 
vuy % — Pazz Y%s) 
+ ag Bay a — 2 Day Ya + Dom dı 
i 2 Dyay ls — Dazu %s + Pyyz %s) 
EN 
+ |, er en ) + u L er +5 WyY 


—, Ws — y) + 4%, %y, + wy(y — %)] 
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Man erhält so a a : ‚ ee ' | & 
n)s — rd), — ; a D), : \ j 
2 er % Se a YtW%)+ Ar % Yl—WwzYs + Wy 2 Yl—W nt “ 


(Us %ı + ds 


| Ba “ W231) -- Ne ways] 


—= ! n « (Bu) —#Ps)s FM. 1a D,— un (War %5 + RR: %g Ys + “yv Y) ee 


—uw + 5 (wit wi) % 


An (BB, Pine z (A DD), — Wn (Ws + #n) RE #.(w} — w}) 5 


x 


Also schließlich 1 
ee (Orden ae ee ee ) 


Damit ist gezeigt: vier Funktionen «, v, ©, w, die den Gl. (2) genügen, erfüllen auf dem 
Rande die Gl. (3), (4), (5). 

3. Wir wollen nun umgekehrt Zeigen: es seien w(x, Y) und D(x,y) zwei Funktionen, 
welche auf dem Rande denjenigen Gleichungen genügen, .die aus (3), (4), (5) entstehen, in- 
dem man auf der linken Seite v=%&, v—= Ö einsetzt, unter = %4(s), % = Ö(s) gegebene 
Funktionen von .s verstanden. Zu w(x,y), D(x,y) gebe es u(#,Y), v(x,y), so daß diese 
vier Funktionen ein Lösungssystem von (2) bilden. (Wenn ©, w der ersten der beiden GI. (1) 
genügen, dann gibt es immer solche « (x, y), v(z, y)). Behauptet wird: es gibt drei Konstante 
Co €, €, derart, daß erstens die vier Funktionen u(2,Yy) + 004 + 6, v(%,Y) — 6% + Cy 
D(x,y), w(xz,y) wieder dem System (2) genügen und daß zweitens längs des Randes -v(x, y) 
+0y+%.=% vlYy)— 69% +0; = © wird. Der erste Teil der Behauptung ist selbst- 
verständlich. Um den zweiten Teil einzusehen, bedenke man, daß die Gl. (3), (4), (5) sowohl 
gelten, wenn a a der linken Seite w, v als auch, wenn man dort %, ö einsetzt. Nach (5) 
also ((ü s—U)m + (% — %)Yyn); =0. Da nach (4) der Ausdruck (%, — wu) + (&, — v,) Yn 
auch in den Ecken stetig ist, gibt es eine Konstante „sodaß (, — u) + (5 Uu)Mm 6% ı 
längs des ganzen Randes. Wegen (8) wird (2, — s)%s = (9) y,=0, also a — u, = ky, 
.- u, = kx, wobei k= c, sein muß. Daraus 4 —u= 00% + co, ds -Vv=— Hr H+% 
wie behauptet. 

Die Gleichungen, die wir für ®, w hergeleitet Haben, insbesondere (3), (4), (5) und die 
ursprünglichen Gl. (1) bleiben alle bestehen, wenn man zu ® eine beliebige lineare Funktion 
hinzufügt. Deshalb darf man für ® zusätzlich fordern 


[J9-42.y = 0 I/®, dedy=0 [Jo azay=0_ RE) 


4. Zusammenfassend können wir sagen: Bei gegebenen RE Nadia sind die 
von Kärmänschen Gl. (1) durch die beiden folgenden a von Randbedingungen 
zu ergänzen: 

Erste Gruppe, Randbedingungen für w: Slider v_(0, w=0 
(eingespannte Platte) oder w= 0, Aw— (1 —») (ws + x%wn) =0 (momentfrei gelagerte 
Platte) oder sonst ein Satz von Randbedingungen, wie er aus der Plattentheorie geläufig ist. 

Zweite Gruppe, N ie Es für DB und w: 


en 


((Dr); IE xD,); ee; 7 (4 Ad. = = WnWg + < 5 (wi + w;) + (u, = DR Yn)s Er (7) 
r +v 
E 


„40 But) ut ty EEE RR EL - (8) 


und in den Ecken s=8;: 
ARM s=,;+0 1 
E (( (Ds —#®,), a 


Außer diesen beiden Gruppen von Randbedingungen verlangen wir für ® zur Normierung 
noch die Gl. (6). 


rl aa 


® $ 
Setzt man En? (h Dicke der Platte), so genügt = @(x, y) folgenden Gleichungen. 


2 - % N 
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Im Innern: - 


N) N Eh® 
AA PR N D p(%, 9 7 7, (Way — Wer Wyy) BD N= 12 (1 2) D 
Auf dem Rande: | | 
Ä | P(s RE BER. 
(An (4) (ra, Pi= 7 (nit tet), 
M (s N /1 
ze) (Pos + % Pn) = ; u)=7(zuitnatun). 
: s=s;+0 P; ei "N/1 s=3,+0 
FÜRTH B=;\a BETEN 127 PER 
schließlich den Normierungsgleichungen: Re 
mdedy=0 Smaray-—0 — ffpdedy=0. 
ER Ä ur K. 


Man kann also das Gleichnis aussprechen: = _ ist die Durchbiegung einer fingier- 


ten Platte, die der Querbelastung p (x, y), der Randkraft P(s), dem Randmoment M (s) 
und den Eckenkräften P, ausgesetzt ist (und die im unbelasteten Zustand mit der gegebenen 


Platte geometrisch zusammenfällt). Elastizitätsmodul Z und Dicke h sind für beide Platten _ 


gleich, hingegen ist die Poissonsche Konstante der fingierten Platte gleich dem Nega- 
tiven der Poissonschen Konstanten der gegebenen Platte®). Aus diesem Gleichnis er- 
gibt sich, daß 9 und damit auch ® bei gegebenem w(x,y) durch die gestellten Bedingungen 


eindeutig bestimmt ist. In der Tat benötigt der Eindeutigkeitssatz der Plattentheorie 


nicht, daß die Poissonsche Konstante positiv sei, es genügt, daß sie zwischen —1 und 
+1 liegt; das ist aber bei der Gleichnisplatte der Fall. 

‚Die Randbedingungen für die linearisierte Theorie ergeben sich, indem man die (in w) 
quadratischen Glieder streicht. Es wird dann @ eine Bipotentialfunktion und in die Rand- 
bedingungen für geht w nicht mehr ein. Das Gleichnis behält auch für die linearisierte Theorie 
seine Bedeutung und ist in diesem Falle ein Gleichnis für die Air ysche Spannungsfunktion 
des ebenen Spannungszustandes bei vorgeschriebenen Randverschiebungen. 

In zwei praktisch wichtigen Fällen heben sich auch in der nichtlinearisierten Theorie 
die in w quadratischen Glieder aus den Randbedingungen heraus. Nämlich bei der einge- 
spannten Platte mit beliebigem Rand (w= 0, w,=0) und bei der freigelagerten Platte mit 


geradlinigem Rand (w=0,#*=0). In diesen beiden Sonderfällen sind also wohl die beiden 
‘Differentialgleichungen (1) für Ö und w nichtlinear, die zugehörigen Randbedingungen aber 


sind linear. 


Die Grundgleichungen für- elastische Platten veränderlicher 
Dicke und großer Ausbiegung. 
Von Karl Federhofer in Graz. 


Die erstmals von Th. v. Kärmän!) entwickelten Grundgleichungen für elastische Platten 
mit großer Ausbiegung gelten für gleichbleibende Plattendicke; es sind dies bei beliebigem 
Plattenumriß zwei partielle nichtlineare Differentialgleichungen 4. Ordnung, die bisher nur im 
Sonderfalle einer am Rand£& eingeklemmten Kreisplatte von St. Way?) und für die am Rande 
drehbar gelagerte Kreisplatte vom Verfasser?) mit Hilfe von Reihenentwicklungen integriert 
worden sind. Im folgenden sollen diese Grundgleichungen unter der Annahme einer beliebig 
veränderlichen Plattendicke aufgestellt werden. Die dabei gegenüber dem Falle gleich- 
bleibender Plattendicke entstehenden zusätzlichen Teile in beiden Gleichungen weisen eine 
bemerkenswerte formale Übereinstimmung auf. 

1. Die Grundgleichungen in rechtwinkligen Koordinaten. Bezeichnen u, v die 
Verschiebungen eines Punktes der Mittelfläche der Platte in den Richtungen eines in der an- 


+) Alle diese Behauptungen liest man direkt ab von den angegebenen Randbedingungen für 9, indem man sie Ver- 


; gleicht mit den Randbedingungen ‚der gewöhnlichen Platte, vgl. z.B. Fr. A. Willers: Die erste Variation der Form- 


änderungsarbeit ausgebeulter ebener Platten. Z. angew. Math. Mech. Bd. 20, S. 118—121, insbesondere Formel Ia und 
IIa, S. 120. 


1) Enz, d. Math. Wiss., Bd.4 (1907—1914), Heft 27, 8. 350. 
2) St. Way: Trans. Amer. Soc. mech. Engr., Bd.56 (1934), S. 627. 
®) K. Federhofer: Luftf.-Forschg., Bd.21 (1944), S.1. 
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} : y Se i 
fänglich ebenen Mittelfläche gewählten «y-System und w(zy) die Durchbiegung der. Platte 
an dieser Stelle, so erfährt ein Element dz dy der Mittelebene die Dehhungen 


a al hu; airs 
= det 3 dx Pk ink SEN ; . Ki RS (1) 
ee 
yTay  2\0y Buch FR “N 
und eine Winkeländerung RR! 
| at ou 00V dwow. u Dr 
Yan 3 9a de y Sk her Sa 0) 9 Be (3), 


wobei die von der großen Ausbiegung w herrührenden Anteile zweiter Ordnung berücksichtigt 

‘ sind, während die in w und » quadratischen Glieder vernachlässigt wurden. Für die diesen 
Größen entsprechenden, über die Plattendicke gleichmäßig verteilten-Spannungen in der Mittel- 
fläche, die wir zweckmäßig zu den auf die Dicke h der Platte bezogenen Spannungskräften 

T,=ho, T,=ho,, S=hr,, zusammenfassen und kurz als Membrankräfte oder — weil sie 
durch die Wölbung der Platte hervorgerufen werden — als Gewölbespannungen bezeichnen 
wollen, gilt nach dem Hookeschen Gesetze‘ » j 


2 | =] “tas) rare nr RN MER . (4) 
ee ea OR a Me Fa SWEET LIEB (5) 
ya 4 yYy ENTE FIRE BUT NEE 0 RN 


. mit u als Querzahl des Plattenstoffes, E und @ als Elastizitätsmodul für Zug und für Schub, 
wobei E=2G11-+u). . BEE 

Da nach Voraussetzung die Durchbiegungen w gegen x und y kleine Größen sind, so 
gelten für die genannten Spannungsresultanten die Gleichgewichtsbedingungen des ebenen 
Spannungszustandes . 


PET oT, 238 ET E / 
| tag ee BE hr RR (7) 
‘ die durch Einführung einer Spannungsfunktion F(xy) gemäß 
ehe ee es ans | 
x y2’ u 9y2® Amy a rn ER = \* (8) 


erfüllt werden. Um die Differentialgleichung, der F zu genügen hat, zu finden, beseitigen wir, _ 
die Verschiebungen u, v aus den Gl. (1-3), indem wir bilden Ri 5 


028, 02 £, | 020 ) 0? w FERN 
ayt Dean ör \noyl Da re (9) 


und hierin die &, &, y,, durch die Spannungsresultanten mit Hilfe der aus (4—6) folgenden 
Beziehungen en 


1 
BEREIT S ae (10) 
| 
Em es (11) 
& » 
/ 1 
Yy= an N EL Er a (12) 


ausdrücken. Dabei wollen wir die Plattendicke A als mit z und y veränderlich voraussetzen. 
Das erfordert zwar eine etwas langwierige Rechnung, die aber zu einem recht einfachen Ergebnis 


führt. Setzen wir vorübergehend 1/R=g, so liefert di k 
und (8) zunächst Ih=g, so liefert. die ei (9) bei Beachtung der Gl. (10—12) 


1 02 92 
Ss 2_ —— —— — 0? 
9 at a For u Fler 
09 9 | ag d 
BE NIE ER 09 8 
og ve BER) H2, +2 re, Fer 


02g N 8) 
tz a u ta ua) +2) = R 
Y 


0 0 
schen Operators A= 7 Su 
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worin die beigefügten Zeiger Ableitungen nach x oder y andeuten. Mit Einführung des Laplace- 
2 2 ı | 
und Ausnutzung der Beziehung 


90 
PE7 ni 


oy° Bw 
AAm)=gAAr + AgAr + 2(8 ar 4302 ö ar) 


gelingt es nach einigen Umformungen, die vorstehende Gleichung 1% die einfache zu 


bringen 


1 1 vd BR | Ber, 
a(zar)+a +mo( 2) u) rer) - 
womit die erste Grundgleichung in allgemeinster Form entwickelt ist, Für konstante Plätten- 
dicke A vereinfacht sie sich in die an bekannte, erstmals von A. Föppl?) abgeleitete erste 
Grundgleichung 
AAF= Eh (w,y? Eu) ET 


Als zweite Grundgleichung für die Funktionen F und tritt hinzu die Biegungsgleichung der 
in ihrer Mittelebene gespannten und durch eine Querbelastung p (je Flächeneinheit) ‘belasteten 
Platte, die bei veränderlicher Dicke lautet 


A(DAw) + (1— pm) [2 Dyywey— Dyylae — Dosty) = 


'f9 ‘9 0 0) 
ie = (Fyu We) Ey 2y (Fz2%y) EP (Fz,Wwy) Bar Fe) EERSFIS (14), 
\ Eh? 2 

wobei D= Di) die veränderliche Plattensteifigkeit bedeutet und der in eckige Klammern 
gesetzte Ausdruck rechts vom Gleichheitszeichen jene zusätzliche Belastung der Platte angibt, 
die von den Spannungsresultanten in der Mittelfläche infolge der durch die Plattenbiegung ver- 
ursachten Schiefstellung der Plattenelemente gegen die ursprünglich ebene, Mittelfläche erzeugt 
wird). Der links vom Gleichheitszeichen in eckigen Klammern stehende Ausdruck läßt er- . 
kennen, daß bei seiner Herleitung die Plattensteifigkeit D als veränderlich vorausgesetzt wurde®). 
Der in eckigen Klammern rechts vom Gleichheitszeichen stehende Ausdruck vereinfacht sich 
nach Ausführung der Differentiationen und Beachtung der Gleichgewichtsbedingungen (7) 
zu P,3Wyy+ FyyWn — 2Fry Wr, so daß die Grundgleichung (14) im Falle konstanter Steifigkeit 
übergeht in?) ’ 
DAAv—=p-+ (FarWyy —2Fay Way + FW): : > (14a) 


Bemerkenswert ist in den allgemeinen Grundgleichungen (13) und (14) die nahezu vollkommene 
formale Übereinstimmung der linken: Gleichungsseiten, die ineinander übergehen, wenn 1/h 
durch D, F durch w und +u durch —u ersetzt wird. 

Mit acer durch Integration der beiden allgemeinen Grundgleichungen gewonnenen Kenntnis 
von F und w sind die Formänderungen und Spannungen in Solchen Platten mit großer Aus-_ 
biegung bestimmt. Die dem Biegungszustand entsprechenden, mit dem Abstande von der 
Mittelfläche linear anwachsenden ERSNIDASN, lassen sich in bekannter Weise zusammenfassen 
zu den Biegungsmomenten 


er = — D (w,r - 4Wy,) RE HNN - a ae = (15) 
d,=—D(wytUW) - : Kr (16) 
bzw. zu einem Scherungsmoment | 
Gu=—-(1- Dun... rennen ae)! 
nnd zu den Querkäften 
} 
N,=,,[P4u) — (1—u) Du, Dive 2 (18) 
0 
MER =, [DAwW— (l—u) ID, w.—D, wu) 2: -: 77 (19) 


"Als Grenzfälle enthalten die beiden Grundgleichungen (13) und (14) mit F=0 den Fall der 


biegungssteifen Platte mit kleiner Durchbiegung (Kirchhoffscher Fall) und mit D=0 den 


Fall der vollkommen biegsamen Platte (Membran). 


“) A. Föppl: Vorlesungen über technische Mechanik, Bd.V. Verlag Teubner (1907), 8.139, Gl. 165. 

5) Vgl.hierzul.B. BiezenoundR. Grammol, Technische Dynamik, 8.584, Gl. (7), Berlin 1939, u.K. Feder- 
hofer u. H. Egger, Ing.-Arch. 14 (1943), S.155. 

°) Siehe H. Marcus: Theorie elastischer Gewebe, Bd.I, S.94, 2. Aufl. Berlin 1932. 

?). In dieser Form, zuerst angegeben von Th.v. Kärmän in Enz. d. Math. Wiss., Bd. 4 (1907—1914), H. 27, S. 350. 


Ir 
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2. Die Grundgleichungenin Polarkoordinaten. Da sich bei den Kreis- und Kreis- ei 

ringplatten, die für die praktischen Anwendungen eine sehr große Bedeutung haben, die Ver- : 
wendung von Polarkoordinaten (r,9) empfiehlt, so wollen wir die Grundgleichungen (13) und (14) 
auf Polarkoordinaten umschreiben. Denken wir uns die X-Achse in die Richtung r, die‘ Y-Achse 

.in die dazu senkrechte Richtung s gelegt, wobei ds= rdg ist, und beachten, daß dann 

| BE SO R ER 

” 98 ra ror’ Orös Or\r 29). 

wird, so lassen sich gemäß (8) die radiale und tangentiale Normalkraft 7, und 7, sowie "die 

Schubkraft 8 durch eine Spannungsfunktion F(r,p) in folgender Form darstellen fe 

i 1.@R7 1108 02F ° 2 

Be 0 — —E an En ee ra ya 20). 
‚rag r dr. or’ 5 or\r dp ER 

. Hiermit lauten dann die Grundgleichungen (13).und (14), wenn der Kehrwert 1/h der Platten- 
dicke zur Abkürzung‘ mit g bezeichnet wird, in Polarkoordinaten 


BEN \[o 2 (1 2g\ 2 Lee ie 4.2 
au 4+ tm), galt rd 


iR 


r op 
ee anl-n N > a 
ar? nee Bi 5 r 99 dr: \r20gp? r or 
o/10D\o /1ow enge on 
DE AH Sr rTT | 
®w/1®2D 10D A ow 0) gw 
Eee rn 2) ee) . . (14b). 


: 0 ow 0 ow 

en or (8 Re (8 =) 
Diese Gleichungen sind erheblich verwickelter als die ihnen bei der Benutzung rechtwinkliger 
‚Koordinaten entsprechenden Gl. (13) und (14), sie erfahren aber eine wesentliche Vereinfachung 
in dem für die praktischen Anwendungen besonders wichtigen Falle einer rotationssym- 
metrischen Belastung und einer rotationssymmetrischen Veränderlichkeit der Dicke h; dann 
wird die Spannungsfläche zu einer Umdrehungsfläche, d.h. # wird unabhängig von 9,'das 
gleiche gilt für die Größen w,g=1/h und D und es gelten unter diesen Voraussetzungen die 
vereinfachten Grundgleichungen ’ 


1 1drd /1\dF E1l dd /dw\? 
A en Dekan ee WERE (136) 
1 d [’@Ddw ; 1d /dFd 
nt) ee ee (14c) 


Beachtet man, daß 


DA o) 


1dıd 
Ada = gr; 


und analog | 
1 1d d /1 
az An) lair)| 


ist, so folgt aus vorstehenden Gleichungen nach einmaliger Integration 


alt ) Ae E (dw\? 

ern a On N. a, 
d dDdw dF dw | | 

ame prar + in. c, Wise nel a (14d). 


Um die mechanische Bedeutung der Integrationskonstanten C', festzustellen, beachten wir, 
daß im vorausgesetzten Falle der rotationssymmetrischen Belastung und ebensolcher Ver- 
änderlichkeit der Plattendicke A die in der Meridiannormalen wirkende Querkraft N gemäß 


den Gl. (18) und (19) gegeben ist durch \ 
d dD1dw 
N=-—[DAw — (1=-u) — — — 
Fe a dr r dr’ 


de j 1dF 
womit sich die Gl. (l14d) wegen 7, = lea umformen läßt in r (nr, nn = [prar +0; 
r 


u 


a A 


ee 3 Are a a rt je 
a nn ” er IE RU 
z en ’ er 2 
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diese Gleichung drückt für die Richtung der Rotationsachse das Gleichgewicht der Kräfte aus, 
die an einem aus der Platte’ herausgeschnittenen Ringe ‚angreifen. | 

Bei gleichförmigem Drucke p ist daher im Falle einer Vollscheibe die Konstante Ur 
gleich Null, während sie z. B. bei der Rin gplatte mit kräftefreiem Innenrand r=r; den Wet 
= erhält. 

Die Integrationskonstante G% ist im Falle der Rotationssymmetrie gleich Null zu en: 
wie die folgende unmittelbare Herleitung der Gl. (13d) zeigt. 

Die elastische Fläche ist eine Umdrehungsfläche, deren Verzerrung durch die radiale 


Verschiebung v(r) und die Neigung p(r) der Meridiantangente gegeben ist, wobei ez = ‚tgp ist. 


Dann gilt für die Ben gen &, und & mit sin ‚us =u 


dv v. 
= dr a EN = r? In 
woraus sich zufolge T, = Fur (&, + u8%) T,= I (& + ue,) die Werte , 
‚vv 1% —uT dv u? 1.7, T, | BR i 
vn Feen und re ern ergeben; hiermit .folgt 


d er E 
A ee | N a a. 


Diese Gleichung läßt sich aber mit Ber der Gleichgewichtsbedingung 


eye Ta=0 
Be dlı \ ) an E 
f ar a es Br 
umformen in Zur TH T-UHrMTT., 2 \% 5 W. ) 
R 1dF d2F | 
Nun ist > SR Er 


womit sich unmittelbar ergibt 


FE RE An ee} 


übereinstimmend mit Gl. (13d), wenn dort C,=0 gesetzt wird. 


Eine Verallgemeinerung der Formel von Pollaczek-Geiringer 
und ihre Anwendung auf- die Verteilungsfunktion 

der Hagelschäden. | 

Von Rolf Lüders in Berlin. l 


Die Formeln der Wahrscheinlichkeitsansteckung bei kleinen Wahrscheinlichkeiten, ins- 
besondere die Formeln von Eggenberger und Pollaczek-Geiringer, sihd 
mit Erfolg zur Berechnung von statistischen Verteilungsfunktionen benutzt worden!). So 
konnte die Verteilungskurve der Hagelschäden unter Verwendung dieser Formeln berechnet 
werden.‘ Dabei wurde angenommen, daß die Schadenswahrscheinlichkeiten konstant sind 
und sich von Jahr zu Jahr nicht ändern. Wie aber die Statistik der Hagelschäden zeigt ist 
‚diese Annahme nicht zutreffend. ‘Es wird nämlich ein periodischer Schadensverlauf beob- 
achtet, wenn auch .infolge der starken Zufallsschwankungen die Periodizität nicht immer 
deutlich zu erkennen ist. Die bekannten Methoden der Periodenanalyse lassen eine Periode 
von 10 bis 11 Jahren entsprechend der Sonnenfleckenperiode erkennen?). 

Im folgenden soll nun eine Formel hergeleitet werden, die die durch die Periodizität 
verursachte Inhomogenität im Schadensverlauf berücksichtigt. An einem Beispiel soll dann 
gezeigt werden, wie sich eine Periodizität auf Grund der Periodogramm-Methode feststellen 
läßt und sich auf die Verteilungskurve der Schäden auswirkt. 

ı) Riebesell: Einführung in die Sachversicherungs-Mathematik. Berlin 1936, S. 48. 


2) Vgl. Myslivec: Periodische' Schwankungen der Schadensquotienten in der Hagelversicherung. Oesterreichische 
Revue, 60.Jg. (1935), S. 379 u.389; Beitrag zur Theorie der Hagelversicherung. Blätter für Versicherungsmathematik, 


3. Bd. (1936), S. 424. 
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‚Zur Berechnung der Wahrscheinlichkeit dafür, daß eine bestimmte Anzahl von Er- 


eignissen in der Zeiteinheit auftritt, dient bekanntlich die Po isson sche Formel, wenn die. 


mittlere Anzahl der Ereignisse je Zeiteinheit nur gering ist. 


Nach dieser Formel ist die Wahrscheinlichkeit dafür, daß genau r Ereignisse.in der Zeit- 


einheit auftreten, 


Bet Brr ae Eee & det nd RT er Bee (1), 


wo m die durchschnittliche Anzahl der Ereignisse in der Zeiteinheit ist. Diese Formel läßt 
sich z. B. auf die Statistik von Krankheiten anwenden. Ist m die durchschnittliche Anzahl 


von‘ Fällen einer bestimmten Krankheit je Monat, so ist P, die Wahrscheinlichkeit dafür, 
daß in einem Monat genau r Krankheitsfälle beobachtet werden. Voraussetzung für die Gül- 
tigkeit der Poissonschen Formel ist aber, daß die Wahrscheinlichkeiten für das Auf- 
. treten der Ereignisse voneinander unabhängig sind, daß also keine Wahrscheinlichkeits- 
ansteckung besteht. : 

Soll auch die Wahrscheinlichkeitsansteckung berücksichtigt werden, so ist die Formel 


von Pollaczek-Geiringer anzuwenden?)*). Ein Spezialfall dieser Formel ist be-. 
reits 1923 von Eggenberger aus dem Urnenschema abgeleitet und auf die Statistik , 


der ansteckenden Krankheiten angewandt worden®). Bei der Formel von Pollaczek- 


Geiringer wird angenommen, daß die Ereignisse entweder einzeln oder zu je zweien oder. 


zu je dreien usw. auftreten können. Ferner wird angenommen, daß die Wahrscheinlichkeit 
‘ dafür, daß ein einzelnes Ereignis bzw. ein Paar von Ereignissen usw. auftritt, von Zeiteinheit 
zu Zeiteinheit konstant ist, daß also Homogenität von Serie zu Serie besteht. Bezeichnet 
man die mittlere Anzahl je Zeiteinheit der einzeln auftretenden Ereignisse mit %,, die-der zu 
je zweien auftretenden Ereignisse mit 2h, und allgemein die mittlere Anzahl der zu je n auf- 
tretenden Ereignisse mit nA,, so ist die Wahrscheinlichkeit dafür, daß in der Zeiteinheit ge- 
nau r Ereignisse beobachtet werden, 


h 4 ( h N) | ar) 
Bi a. I )\.(e-h » 22). Joh.» 
R P,= 2, (e nt i v.| x Y,! 


+21, +3%,4....=r 
[4 , 
P,=eh-hh—.... a 
ahardet.. Sn hi eely ne i 
Die mittlere Anzahl der Ereignisse in. der Zeiteinheit ist 


n=hr2ht amt. 


Für die Praxis ist die obige Formel nicht ohne weiteres brauchbar, weil die A,, ha; has...» 
meistens nicht unmittelbar beobachtet werden können sondern aus der Häufigkeitsverteilung 
berechnet werden müssen. Das ist aber nicht möglich, weil es sich hier um unendlich viele 
Parameter handelt. Es ist daher zweckmäßig, die Formel zu spezialisieren, etwa dadurch, daß 
man die Ah; in einer bestimmten Progression abnehmen läßt, z.B. 


.2). 


2 N—2 
h=5 =, = Fr mn... EEE LS 

Nun ist aber in vielen Fällen die Wahrscheinlichkeit für das Auftreten eines Ereig- 
nisses nicht zeitlich konstant sondern ist mehr oder weniger starken Schwankungen unter- 
worfen. Bei den meisten Krankheiten schwankt z.B. die mittlere Zahl der Krankheitsfälle 
mit der Jahreszeit; ferner wird vielfach bei Beobachtung eines längeren Zeitraums eine Ab- 
nahme der Wahrscheinlichkeit im Laufe der Jahre festgestellt. Die Formel von Pollaczek- 
Geiringer ist nun so zu verallgemeinern, daß sie auch auf die Statistik von Ereignissen mit zeit- 
lich veränderlichen Wahrscheinlichkeiten (Inhomogenität von Serie zu Serie) angewandt 
werden kann. Zu diesem Zwecke sei angenommen, daß ein Zeitabschnitt in n Teilabschnitte 
zerlegt werden kann und in dem &-ten Teilabschnitt ähnlich wie in Formel (2) die mittlere 
Anzahl der einzeln auftretenden Ereignisse A®, die mittlere Anzahl der zu je zweien auf- 
tretenden Ereignisse 2h® und allgemein die mittlere Anzahl der zu je n auftretenden Er- 
eignisse :n : h® ist.‘ Dann ist nach GI. (2) die Wahrscheinlichkeit dafür, daß in dem k-ten 
Teilabschnitt ‚genau r Ereignisse beobachtet werden, 


h%) [21 h\%) v, © 


vl wel 


(k) 


RE N 
r 


ni 


(4). 
Me ee NE ,+29%+...=r 
®) Pollaczek-Geiringer: Z. angew. Math. Mech. 8 (1928), S. 292. 
‘) Lüders: Die Statistik der seltenen Ereignisse. Biometrika 26 (1934), 8. 108. 


5) Eggenberger und Polya: Über die Statistik verketteter Vorgänge. Z.angew. Math. Mech. 3 (1928), S, 279; 
Eggenberger: Die Wahrscheinlichkeitsansteckung. Bern 1924, 
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Für einen längeren Zeitäabschnitt ist dann die relative. TRulblekeir dr Teilabschnitte mit genau ' 
r Ereignissen 


1 
BA de (5) 
k=0 
Die mittlere Anzahl der Ereignisse je Zeiteinheit ist dabei 
n—1 M 
m—=— > m) 2 
ın 
4 'k= 0 N 

wo m® — hd + 2m + 3.1@ N RL Ba 
der entsprechende Mittelwert für den k-ten Teilabschnitt ist. S 


Die Formeln (4) und (5) lassen sich im Prinzip ebenso anwenden wie die Formel (2), die 
für den zeitlich homogenen Ereignisverlauf aufgestellt war. Es möge sich z. B. darum handeln, 
die Verteilungskurve für die monatlichen Erkrankungsfälle einer bestimmten ansteckenden 
Krankheit, für die ein Beobachtungszeitraum vön 10 Jahren vorliegt, zu bestimmen. Dabei 
sollen starke jahreszeitliche Schwankungen auftreten, während von Jahr zu Jahr die Wahr- 
scheinlichkeiten sich nicht ändern. Man hätte dann aus der beobachteten Häufigkeitsver- 
teilung der Januar-Krankheitsfälle die 7%, aus der Februar- -Verteilung die A“ usw., aus der 
Dezember-Verteilung die „ıD zu bestimmen und hieraus die pP (k:==.0, 1, 2,2 u 1. und 
P, zu berechnen. Diese Methode führt aber meistens nicht zum Ziel, da die Anzahl der 
beobachteten Werte (in dem erwähnten Beispiel für jeden Monat nur 10) zu gering ist, um aus 
der beobachteten Verteilung die Parameter der theoretischen Verteilung mit Bendgendet Ge- 
nauigkeit zu bestimmen. ü 

Es muß also ein Weg gefunden werden, der die Berechnung der Parameter aus dem 
Gesamtmaterial, bei dem sich die zufälligen Schwankungen besser ausgleichen, ermöglicht. 


Dabei müssen über die A(®, h®, .... bestimmte Voraussetzungen gemacht werden. Es 
sei daher angenommen, daß die A®, A®, ... dem Mittelwert m(% proportional sind: 
BEN RI BG ne ee SE 
m%) ; 
wobei 0; = = also ml — 0, m ist. 
Die Aufgabe besteht jetzt darin, aus der beobachteten Verteilung die 715 27%, . . . zu bestimmen. 
Dabei sollen entsprechend den Gl. (3) wieder die folgenden Annahmen gemacht werden: 
| Base 
m=5 mL, et ee: a Te BEER A dt (7a). 


Man könnte nun die n, ß, y mittels der ersten drei Momente berechnen. Da aber das 
Moment zweiter Ordnung und insbesondere das Moment dritter Ordnung mit großen mitt- 
leren Fehlern behaftet sind, soll auch hier die Methode der modifizierten Mittelwerte ange- 
wandt werden [vgl. %)]. Zu diesem Zwecke sei die folgende erzeugende Funktion eingeführt 


n—1 


fis;2) - Ir, PR 22 Be N eg eraild 


Dabei ist 
x © hk)”: rk”: ERS 
WM —= Pi) ar Mn... 1 ar zr 
: 2 N v!vz! un 
= ee -M—M—.. Be NETZ az... (9). 
Hieraus erhält man unter Benutzung von (7) 
TE a ua ie 2.00, EN (10). 
Es ist mithin : 
—1 y & 
an) BL n mn. .+maz+maa+..}C, BT (11) 
; rt, 
und 
et # 2% 1 
fix, )=— Yeti" + mit NMm® k = E (x)°% (12) 
n 9 k=0 
wo 
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ist. Differenziert man die Funktion ’f (z, z) nach z und setzt z= 1, so erhält man den Ausdruck 

al) = j@, DM 2 ER 3 (14). 
Andererseits erhält man aus (11) 


n—17 


f(@1)= 2 matt et.) 


En y 
= = (a An I 0, E(z)°%; 


es ist also 
f 1 A Es 
tm rat. er de) er (15), 
— + DO, B(a)'% 
k=0 


wo die rechte Seite der Gleichung zur Abkürzung mit &(x) bezeichnet worden ist. 
Aus den Gl. (7a) folgt 
, x 
mM+2n2+3n +... =m+Pßr+PByr+ By? BEE A MT er 
Nimmt man drei verschiedene Werte für x, nämlich x,, x, und x, an, so ergeben sich die Glei- 
chungen 


m+ „a Zn ale a Te RER (16) 
N ar N (17) 
m+ By ER RS REIF LUREBRS SOE SONS IR SSR, (18) 


Hieraus folgt » | 
1—y2 - 2% —% , dla) —älx;) 


ee (19) 
Bee ya) en nr tie ir. 
ma) HE NTRREN A REN. 


Aus diesen drei Gleichungen können die y, ß, n leicht berechnet werden, falls die äa), ä&(%s) 
und ä(z;) bekannt sind. Nun gilt für &(x) — vgl. (15) — 


Ri 1 2ft%1) 


dabei ist 
el- Drha 
aus den beobachteten Werten der P, leicht zu berechnen. Für E(«) gilt die Gl. (12) 


n—1 


DE rn HER N en a hen. (23), 


"wo f ft = I Bw 


ebenfalls aus den beobachteten P, leicht zu nen ist. Der Ausdruck E (x) ist eine Wurzel 
der algebraischen Gl: (23) und 1äßt sich mittels eines Näherungsverfahrens ‚berechnen. Da 
die linke Seite dieser Gleichung 


| Ei)" + E(a)"+....+ E(a)-1 
unter der Voraussetzung 0,>0 eine monoton wachsende Funktion von E(x) darstellt und 


für Zi) = 0 verschwindet, dagegen für E(x) > » über alle Grenzen wächst, da ferner stets 
/®1)>0 für &>0 gilt, hat die Gl. (23) genau eine reelle Wurzel E(x). 


u Se ae 


\er, 
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" Man’ hat also drei Werte für x anzunehmen, wobei zweckmäßigerweise x, = T’ gesetzt. 
wird, und aus den beobachteten Werten der P, f(z,1) und f’(z,1) zu bestimmen. Aus der 
Gl. (23) erhält man dann E(x,), E(x,), E(z,) und aus der Gl. (22) ä(z,), ä(x), &(x,). Die 
Gl. (19), (20) und (21) liefern die Parameter n,, ß und y, mit deren Hilfe dann Ken Benutzung 
von (7) und (7a) aus den Gl. (4) und (5) die theoretische Häufigkeitsverteilung ermittelt 
werden kann. Hierzu werden nur noch die Werte der CO, benötigt. Die Bestimmung könnte 
aus den beobachteten Mittelwerten m® der Teilabschnitte mit Hilfe der Gleichung 
m * 


erfolgen. In ddm besprochenen Beispiel hätte man also aus den 120 beobachteten Werten den 
Mittelwert m, aus den 10 beobachteten Januar-Werten m, aus den 10 Februar-Werten m) 
usw. zu bestimmen. In der Regel wird dieses Verfahren aber ‘deswegen nicht anwendbar 
sein, weil die Zahl der Werte, aus denen m(%) bestimmt wird, zu gering ist. Deswegen er- 
scheint es zweckmäßig, für den Verlauf der m% eine bestimmte Funktion anzunehmen und 
unter Benutzung dieser Funktion die m®) festzulegen..Bei der Statistik von Krankheitsfällen, 
deren Zahl mit der Jahreszeit stark schwankt, RD es sich, einen sinusförmigen Verlauf 
der m(% anzunehmen: ; 

m) —= m + cos ei bsin en 


=m-+ Asin (yo + = 


wobei die a und 5 bzw. A und %, aus den beobachteten Mittelwerten mittels der Methode 
der kleinsten Quadrate zu berechnen sind. Auch im Falle eines nicht periodischen Verlaufs, 
z. B. eines linearen oder parabolischen Trends, kann man die verallgemeinerte Poisson- 
sche Formel anwenden und die theoretische Verteilungskurve aus den Gl. (4) und (5) bestim- 
men. Man hat die Konstanten der Trendgleichung mittels der-Methode der kleinsten Quadrate 
zu bestimmen und dann die h(® zu berechnen. 


(24), 


‚ 
t 


ee 


Für eine Reihe von Infektionskrankheiten ergibt sich ein dsutlicher periodischer jahres- 
zeitlicher Verlauf der Erkrankungs- bzw. Sterbefälle). Es liegt daher nahe, als Beispiel für 


die Anwendung der oben hergeleiteten Formel die Häufigkeitsstatistik einer Infektionskrank- 


heit zu benutzen. Es sollen die wöchentlichen Todesfälle durch Meningitis (epidemische 
Genickstarre) in der Gesamtheit der deutschen Großstädte in den 10 Jahren 1925 bis 19347) 
untersucht werden. Um die Feststellung des periodischen Verlaufs zu erleichtern, werden 
jeweils die Todesfälle für 4 aufeinander folgende Wochen zusammengefaßt, so daß das Jahr 
in 13 4-Wochen-Abschnitte zerfällt. Der Mittelwert beträgt m’ = rd. 10,6 Todesfälle je 4- 
Wochen-Abschnitt. 


Tafel 1. Der periodische Verlauf der Todesfälle durch Meningitis. 


III. IV. 

II. Beobachteter Theoretischer ' 

I. Mittelwert f. d. Mittelwert f. d. 

k k-ten 4-Wochen- | k-ten 4-Wochen- 

Abschnitt Abschnitt 

1.— 4. Woche 0 11,2 12,07 
5 8. Ar 1 14,2 13,10 
MA 2 17,4 13,55 
13.—16. ss 3 12,4 13,32 
17.—20. $ 4 10,1 12,46 
21.—24. 3 5 11,2 11,16 
28.387"), 6 9,3 9,73 
29.—32. ” 7 9,4 8,49 
33.—36. 57 8 81 7,12 
BT. AD 9 fi! 7,62 
41.—44. ven 10 9,8 8,18 
45.—48, Ar 1l =, 9,29 
Aare 12 7,4 10,70 


°) Vgl.z.B.Ziezold: WeitereErgebnisse der Analyse von Seuchenkurven. Z.Hyg. Infekt.-Krankh., Bd. 124 (1942), 


8.93. 
?) Reichsgesundheitsblatt, Jg. 1—10. 
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Spalte: III den Tafel 1 zeigt die beobachteten Mittelwerte für die einzelnen 4-Wochen- 
Abschnitte; bereits aus: diesen Zahlen ist der .periodische Verlauf deutlich erkennbar. Zum 
Nachweis der Periodizität ist das sogenannte Periodogramm?) zu berechnen. Zu diesem Zwecke 


sind die folgenden Werte für r=3,.r = 4 usw. zu ermitteln: 
. h R ’ 
(gi 


Kurz ı x ER ' X; an 
= > (y, #4 + Mror +... )00s 9 
‚=0 e b \ - 25 
nd (23) 

ER! ERS n 2” 
ra U RE a re 


K} 


Dabei bedeuten N = 130 die Gesamtzahl der 4-Wöchen-Abschnitte und die y, die Ab- 
weichung der im »-ten 4-Wochen-Abschnitt beobachteten Zahl der Todesfälle von dem (ab- 
gerundeten) Mittelwert 10,6. (Die 4-Wochen-Abschnitte sind laufend von 0 bis 129 numeriert. 
Aus den a, und 5, erhält man die Amplitude für eine Periode von r Abschnitten 


A.=Vu+b. 


Es ergibt sich das folgende Periodogramm: 


r | 4, r | 4, r Ar 
3 0,38.:| -8 |; 0,29 13 2,99 
4 0,30 9 1,37 14). 118} 
5 0,40 10. 0,50 15 .| 1,32°_ 
6 033 | ı 0,94 16 1,67 ” 
7 09 | 12 1,09 17 0,42 


Das Periodogramm zeigt ein deutliches Maximum für r = 13, d.h. die Periode 13 ist vor- 
herrschend, womit die behauptete jahreszeitliche: Periodizität erwiesen ist. Man erhält für 


r = 13 a,, = 1,50 und b,, = 2,59; die theoretischen Werte für den k-ten 4-Wochen-Abschnitt 
sind 


2 .. 27 
y=m’ + 1360875 k+ busin nk. 
Die Rechnung führt zu den in Spalte IV der Tafel 1 angegebenen Werten. Diese Zahlen zeigen 
einen ähnlichen Verlauf wie die beobachteten Zahlen, wenn auch die letzteren stärkere Schwan- 
kungen aufweisen. Eine bessere Anpassung würde man erhalten, wenn man auch die höheren 
Glieder der’ harmonischen Entwicklung berücksichtigen würde. Da aber nicht feststeht, ob 
diese Perioden noch reell sind, soll auf die höheren Glieder verzichtet werden. K 


Tafel 2. Die Häufigkeisverteilung der Todesfälle durch Meningitis. 


en 


I. II. II. IV. | V. 
2 Häufigkeit nach Häufigkeit nach der 
Zahl Enter Beobachtete Häufigkeit re 0 Yas BER 
dizität) ohne Berücksichtigung mit Berücksichtigung 
der Periodikität der Periodizität 
a N EFT a FR a a A Er Fr EN 
0 54. 37,0 52,3 54,0 
l 114 97,8 112,3 111,2 
2 120 129,3 124,2 ; 121,9 
3 92 113,9 96,0 95,5 En 
4. 61 75,2 59,8 61,0 
5 32 39,7 331 34,5 
6 25 ITIDugR aber 18,3 
ML 10 6,6 In 97 
8 5 2,2 5,5 5,3 
9 2 0,6 3,3 3,1 
10 2 0,2 2,1 1,9 
11 u. mehr 3*) 0,0 ° 4,0 3,6 
Summe | 520 | 520,0 | 520,0 520,0 


®) Vgl. Stumpff: Grundlagen und Methoden der Periodenforschung. Berlin 1937. 
*) Es wurden je einmal 11, 12 und 19 Todesfälle beobachtet. 
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Tafel2 zeigt nun in Spalte II die beobachtete Häufigkeit (wöchentliche Zahlen), in 
Spalte III die Werte der Poisson schen Formel, die von den beobachteten Werten syste-: 
matisch abweichen, und in Spalte IV die Werte der Formel von Pollaczek-Geiringer, 
die unter Benutzung von lg x, = 0, Ig x, = —0,02 und lg x, —= —0,04 berechnet worden 
sind?) und mit der beobachteten Verteilung befriedigend übereinstimmen. Bei den Zahlen 
der Spalte V ist auch die Periodizität berücksichtigt ROERSE: dabei wurden die a, und b, für 
r = 52 berechnet. 


"= 0,258, bu 0,716 


und hieraus die m® nach Gl. (24) bestimmt. Die Zahlen der Spalte V unterscheiden sich 
nicht wesentlich von den Zahlen der Spalte IV; die pa unR an die beobachtete Verteilung» 
ist eher noch etwas besser. 


Bei der Anwendung der Formel von Pallaczek‘ Gelriefer auf die Hagelstati- 
stik könnten insofern Bedenken auftreten, als die Dinge hier nicht so einfach wie bei der 
Statistik der Unfälle oder der ansteckenden Krankheiten liegen. Wenn auch durch einen 
Hagelschlag ein mehr oder weniger großes Gebiet getroffen werden kann, so daß ähnliche 
Verhältnisse wie bei der Krankheitsstatistik vorliegen, so ist doch die Zahl der während eines _ 
Hageljahres getroffenen Gebiete zu groß, um daraus die Gültigkeit des Poissonschen 
Verteilungsschemas zu folgern. Auch die Abrundung auf volle Promille-Sätze und die Gleich- 
stellung dieser Sätze mit der Ereigniszahl!®) entspricht nicht den tatsächlichen Verhältnissen. 
Andererseits haben aber die durchgerechneten Beispiele gezeigt, daß die erweiterte Poisson- 
sche Verteilungskurve zumindest eine gute Annäherung an die empirische Schadensverteilung 
ergibt, so daß diese Formel vorläufig. eine DESUCOhEen Methode für die Ermittlung der Schadens- 
verteilungskurve liefert. 


Als Beispiel soll die Statistik über die ee Fläche in Württemberg in den 
60 Jahren 1828 bis 1887 untersucht werden"*). 


Tafel 3. Die Häufigkeitsverteilung der Hagelschäden in. Württemberg 1828—1887. 


T: I. HS: IV. V. | VI. 
Verhagelte Fläche Beobachtete Häufigkeit nach a Are iis och De 
im oo Häufigkeit der ZETSspEen Fonmel ea 
B ormel T n Formel von 
Pollaczek- Geiringer Poissonschen Formel. Pollaczek-Geiringer“ 
0— 2 5 0,3 4,8 0,8 4,9 
3— 5 15 5,8 13,4 8,4 N 13,2 
6— 8 14 19,4 14,0 18,0 13,9 
9—11 6 21,3 10,7 17,4 10,7 
12—14 10 10,2 71 10,3 TEL 
15—17 4 2,6 4,3 3,9 4,3 
18—20 2 0,4 2,5 1,0 2,4 
21 u. mehr EBENE REBEL EEE ER USE RE 
Summe | 60 | Eee Fe 


Tafel 3 zeigt in Spalte I die verhagelte Fläche in °,, in Spalte II die Zahl der 
Jahre, in denen die in Spalte I angegebenen Promillesätze beobachtet worden sind. Dabei 
sind die in den einzelnen Jahren festgestellten Schadenssätze auf volle Promille abgerundet 
worden. Mit diesen Zahlen wurde dann die Rechnung durchgeführt; nur der Übersichtlich- 
keit halber sind in der Tafel 3 jeweils die Promillesätze 0—2, 3—5 usw. zusammengefaßt 
worden. Der Mittelwert beträgt m = 9,233, d.h. im Durchschnitt der 60 Jahre ver- 
hagelten jährlich 9,233%/,, der Gesamtfläche. Spalte III zeigt die nach der Poissonschen 
Formel berechneten Zahlen, die. systematisch von den beobachteten Zahlen abweichen. 
Spalte IV zeigt die Werte nach der Formel von Pollaczek-Geiringer ohne Berück- 
sichtigung der Periodizität, wobei z,=1, %, = 0,97724, x, = 0,95500 gesetzt wurde. 
Die Übereinstimmung mit den beobachteten Werten ist befriedigend bis auf den Wert 
für 9—119/- 


») Zur Berechnungsmethode vgl. Lüders: a.a. O. 
„.2) Vgl. Lüders: Über die Verteilung der Schadenquotienten bei Hagelversicherungsunternehmungen. Neumanns 
Zeitschrift, Bd. 57 (1934), S. 761. 

u) Bühler: Die Hagelbeschädigungen in Württemberg während der 60 Jahre 1828—1887. Württembergische 
Jahrbücher für Statistik und Landeskunde, Jg.1883, I. Bd., 3. Heft. 
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Zur Feststellung, ob ‚ein periodischer 'Schadensverlauf vorliegt, ist das Periodogramm 
zu berechnen. Man erhält nach der geschilderten Methode die folgenden Werte: 


x 
u Tr 

3 1,26 12 0,60 f 
4 0,43 13 1,47 
5 1,80 14 1,87 

6 0,46 15 1,04 

7 2,08 16 0,99 

8 0,92 17 1,02 

9 2,74 18 ‚| 1,02 ra 

10 2,86 ..19 1,08 

ir\.40. 01,70 20 | 0,80 


Es ergibt sich ein deutliches Maximum für r = 10, d.h.die Schäden zeigen eine angenähert 
10-jährige Periodizität. Man erhält a, = —2,52 und b,, = 1,36. 

Wie die Spalte V der Tafel 3 zeigt, ergibt die Poissonsche Formel auch dann noch 
keine befriedigende Annäherung ‚an die beobachtete Häufigkeitsverteilung, wenn man die 
Periodizität berücksichtigt. Nach der Formel von Pollaczek- Geiringer „(Spalte VD) er- 
hält man fast dieselben ‚mit der Beobachtung bis auf den Wert für 9—11°/,, befriedigend 
übereinstimmenden Zahlen wie in Spalte IV. 

Bei dem behandelten Beispiel fällt auf, daß die nach der ‚‚unperiodischen“ Formel 
berechneten Werte nur wenig von den Werten der ‚periodischen‘ Formel ‚abweichen. 
Der Grund für diese Erscheinung ist darin zu suchen, daß die P, der „periodischen“ 
Verteilung als Summe von % Gliedern dargestellt werden, die nach der Gl. (4) fberechnet 
sind, wobei man die Parameter aus den 71, N .... durch Multiplikation mit Oz. er- 
hält. Diese Summe ergibt nun bei nicht zu großer Amplitude angenähert denselben Wert wie 
die Anwendung der Formel (2) ohne Berücksichtigung der Periodizität. Die Annäherung 
wird noch dadurch verbessert, daß die Parameter der Formel aus den modifizierten Mittel- 
werten der empirischen Verteilung berechnet werden. 

Wie die Untersuchung von weiterem Material der Hagelschadenstatistik, insbesondere 
der Schadensergebnisse einzelner Versicherungsgesellschaften, zeigt, wird fast immer ein peri- 
odischer Schadensverlauf mit einer Periode von etwa 10 bis 11 Jahren beobachtet. In vielen 
Fällen fallen die Jahre mit den größten Schäden mit den Jahren, in denen ein Maximum 
der Sonnenflecken-Relativzahlen beobachtet wurde, zusammen oder liegen jedenfalls in deren 
Nähe. ’ 

Bei der Berechnung der Rücklage. für den schwankenden Jahresbedarf nach der mathe- 
matischen Methode wird es sich, wenn der periodische Schadensverlauf für die betreffende 
Anstalt nachgewiesen. ist, empfehlen, diese Periodizität zu berücksichtigen. Die Rücklage 
wäre dann in eine Rücklage I und eine Rücklage II aufzuspalten. Die Rücklage I berücksich- 
tigt die Schäden, die sich aus dem periodischen Schadensverlauf ergeben; sie wird aus den 
ersparten Nettoprämien während der Minimum-Jahre gebildet und in den Maximum-Jahren 
wieder aufgelöst. Die Rücklage II wird für diejenigen Überschäden gebildet, die sich durch 
zufälligeSchwankungen ergeben können. Diese Rücklage ist in der Regel erheblich höher 
und läßt sich mit Hilfe der Formel von Pollaczek-Geiringer berechnen!?2), Da die P® 
sich wegen Gl. (4) von Jahr zu Jahr periodisch ändern, würde sich bei genauer Rechnung 
auch für die Rücklage II ein mehr oder weniger hoher Bedarf ergeben, je nachdem ob es sich 
um ein Maximum- oder Minimum-Jahr handelt. Es muß je nach den besonderen Verhält- 
nissen entschieden werden, ob auch diese Schwankungen bei der Ermittelung des Rücklagen- 
bedarfs berücksichtigt werden sollen. 

Die Berechnung der Schadensverteilungskurve für eine einzelne Versicherungsanstalt 
ist nicht immer so einfach durchzuführen wie bei dem untersuchten württembergischen Material. 
Soweit nämlich überhaupt ein längerer Beobachtungszeitraum zur Verfügung steht, ist der 
'Versicherungsbestand häufig nicht homogen geblieben. In einigen Fällen hat sich der Bestand 
z. B. stark vergrößert, so daß die älteren Schadenszahlen nur bedingt verwertbar sind. In 
vielen Fällen wird aber die Bearbeitung des Materials nach der beschriebenen Methode mög- 
lich sein und damit die Beantwortung der Frage, wie hoch eine Rücklage sein muß, damit 


noch Schäden gedeckt werden können, die mit einer bestimmten Wahrscheinlichkeit zu er- 
warten sind. 


‘») Vgl. Lüders: Neumanns Zeitschrift, Ba.57 (1934), 8. 761 ft. 
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Zur Bestimmung der mittleren Abscisse. 


Bei der zeiehnerischen Bestimmung der mittleren 
Abscisse eines Kurvenabschnittes pflegt man so vor- 
zugehen, daß man zur Sehne des Bogens an ihn eine 
parallele Tangente zeichnet, durch den Halbierungs- 
punkt der Sehne eine x-Parallele legt und das Stück 
zwischen Sehne und Tangente in dreigleiche Teile teilt. 
Der der Tangente zuliegende Teilpunkt gibt annähernd, 
die mittlere Abscisse, Verfährt man genau so mit dem 
y-parallelen Stück zwischen Sehnenmitte und Tan- 
gente, erhält man einen Näherungswert für die mittlere 
Ordinate, Wie man leicht sieht, schließen die so be- 
stimmten Stücke der beiden Ersatzstufenkurven in dem 
betreffenden Abschnitt mit der x-Achse flächengleiche 
Stücke ein, y 

Durch diese Konstruktion wird allgemein nicht mehr 
der durch Parabeln dritten Grades bestimmte Flächen- 
inhalt exakt wieder gegeben, wie durch eine Konstruk- 
tion nach der Keplerschen Faßregel!). Der Fehler 
obiger Konstruktion berechnet sich folgendermaßen, 


‚wenn man etwa von der Konstruktion der mittleren 


Ordinate ausgeht. Man legt die Ordinatenachse in die 
Mitte des Intervalles und bezeichnet seine Grenzen mit 
+h und —h; dann hat man zur Bestimmung der Ab- 
scisse ö des Berührungspunktesder zur Sehne parallelen 


Tangente 
PERS t ch) 
3h =f (6) 
oder, wenn man beiderseits um x =0 entwickelt und mit 
den ersten Gliedern abbricht, 


KA+zw+...=nttdt... 
Also ist N 


DZ 


0 
> 
N 0, 


Die Gleichung der Tangente in dem Punkt mit der 
Abseisse ö lautet 
= FNR— 9: 
Also wird die Ordinate der Tangente für x = 0, wenn 
man sich auf die erste Näherung beschränkt, 


et 14 d No 
Y%=F(ö) — 67 (6) hf ö° 
ara 
ef 12f% hr. 
Damit ergibt sich für den Näherungswert 


h 
sywoatiyty4mM 
1 # IL AEH Zu 
=2jit zh + | se Hp|er... 


während der exakte Wert 
h 


+ 

1 (IV) 
[r@az=21%+ Ft _ +... 
—h 


ist. Das erste Fehlerglied ist somit: 


RER a a 
90 ap” [er 
Der Fehler ist auch hier von fünfter Ordnung in A. 
Aber da f’’”’ im Restglied vorkommt, wird eine Parabel 
dritten Grades hier nicht mehr, wie bei der Konstruk- 


'tion nach der Keplerschen Faßregel, exakt integriert. 


Den Fehler der Konstruktion nach dieser Regelschätzt 
das erste Klammerglied ab!). Das zweite gibt den 
durch die angegebene Konstruktion hinzugekommenen 
Fehler an. Dieser wird dann klein, wenn f, groß gegen 
{= ist, d.h. wenn der Bogen nicht allzusehr von einem 


2) Fr. A. Willers: Methoden der praktischen Ana- 
lysis. Berlin 1928, 515 (22). 
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KLEINE MITTEILUNGEN 


gewöhnlichen Parabelbogen abweicht. Keinesfalls ist 


also die Konstruktion erlaubt, wenn die Kurve im In- 
tervall einen Wendepunkt hat. Haben fIV und f’’ das 


gleiche Zeichen, und ist absolut genommen das zweite 


Glied nicht größer als das erste, so wird der Fehler 
kleiner als bei der Keplerschen Faßregel. Genaue Werte 
liefert diese Konstruktion immer dann, wenn es sich 
um. den Bogen einer Parabel handelt, auch wenn ihre 
Achse .der y-Achse nicht parallel ist, also z.B. für 
Kurven, deren Gleichung y = Ax + BYx+ CO lautet. 
Dresden | Fr. A. Willers, 


Die Querschwingungen der Schraubenfeder., 
1. Bezugnahme auf das Schrifttum. Die Arbeit 
J. Dick, The transverse vibrations of a helical spring 
with pinned ends and no axial load, Philos. Mag., 
J. Sci.,, London, 7, 33 (1942), 513—519 soll berichtigt 
werden. wi \ { 
2. Aufstellung der Differentialgleichungfür 
die Querverschiebung. I 
x Bogenlänge der Schraubenachse, 


R 4Q m auf die Längeneinheit der 


bett Schraubenachse bezogene 
Bee er | M Masse, Ku 
I} w Querverschiebung. ; 


Die Bewegungsgleichung: Kraft = Masse mal Be- 
schleunigung ergibt 


0x I te 
D = Durchmesser der Schraubenlinie. 
Das auf die Längeneinheit der Schraubenachse be- 
zogene Massenträgheitsmoment ist nach Dick, S. 514, 
2 


Gl. (1) m = Die Bewegungsgleichung: Kraftmoment 
— Trägheitsmoment mal Winkelbeschleunigung ergibt 
u) 


+95" dem 
v = Verhältnis der Querkürzung zur Längsdehnung, 
F = Längssteifigkeit, 
l = Länge der Schraubenachse. 
Die Schersteifigkeit ist nach Dick, 8. 514: 2 (1+») Fl. 
w, durch die Querkraft verursachte Querverschiebung 


dw, 
TEE En a A a). 
Die Biegungssteifigkeit ist nach Dick, $. 514 
1495, D? ; 
2+v 2 


w, durch das Biegungsmoment verursachte Querver- 
schiebung 
1-+» D: 0? wg \ & 
M=37, age na: 
Die ganze Querverschiebung ist w= w,+ w,. Dick 
setztin der Bewegungsgleichung Kraftmoment = Träg- 
heitsmoment mal Winkelbeschleunigung statt der 
ganzen Querverschiebung nur die durch das Biegungs- 
moment verursachte Querverschiebung ein, obgleich 
er in der Bewegungsgleichung Kraft= Masse mal Be- 
schleunigung die ganze Querverschiebung einsetzt. In 
den Gleichungen, die die Querkraft und das Biegungs- 
moment durch Differentialoperatoren der Teilverschie- 
bungen ausdrücken, werden die Teilverschiebungen 
eliminiert 
1 Cs a _R"w 
3(l+v)Flöx 1+4»FiD: PEN 

In den Bewegungsgleichungen und der letzten Gleichung 
wird die Querkraft eliminiert 
w Mm _„D o®w 
ae de 8 0m0n 

m ww 24» 2 Ar w 
2(1+»)FI dt 14» FiD: da: 


B)- 
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In den letzten beiden Gleichungen wird das Biegungs- 
moment eliminiert u 
ww At m Aw 24» m. Van, 
om Im drideoer "IH FImBom 
3. Integration der Differentialgleichung für 
die Querverschiebung. Der Bernoullische Pro- 
duktansatz Ww=A(z)sin (wo? -+ a) ergibt diegewöhnliche 
_Differentialgleichung ; 
dA 44+v m „dA 2+v m 
det "Ir Ari” de IHrrD 
Ihr allgemeines Integral ist 


aA=0. 


A= K,SinA,x + Kg008A,x+ K,sinA,0+ K,cosA,z, 


worin } 
ET Rue 
OF], 8r® 
| Sr ER N A La E 
% VG Fa ray, FD" 

4. Erfüllung der Randbedingungen. Sind die 
Enden gelenkig gelagert, so ist dort w=0 und M=0. 
Nach GI. (3) wird die Randbedingung M =0 mit Hilfe 

i | 
der’Randbedingung w=0 zu = 0 vereinfacht. Die 
Enden der Schraubenachse werden als Punkte 2=0 
und #=]/ gewählt. Die Bedingungen am Rande 2 =0 


ergeben 
N K,=0 K=%; 
die Bedingungen am Rande «=1 ergeben 
K,=0 Al=nn Mal, 2, 3,12. 
5. Ergebnis. Für die Eigenfrequenzen gilt 
De n‘ ni F 
ll 4) ——— — 
(4-+»)n?n?+8(2-+) Di 
Die ganze Querverschiebung ist 


ml 


ENTER. 
j w= K,sin 7 sin (w!-+.o). 
Nach Gl. (3) ist das Biegungsmoment 
Re 12 ; 
ME a Te 
Ei 3 3 W. 
) d+rn)n’n® +82 +) Dr 
Nach der Bewegungsgleichung: Kraftmoment = Träg- 
heitsmoment mal Winkelbeschleunigung ist die Quer- 
kraft ; 
3 73 Ö 
Q=4(1+) uud Rt. 
’ “rnmm+82 +) 
"7=sin (wt +0). 


Aus Gl.(1) folgt die durch die Querkraft verursachte 
Querverschiebung 
2 n?n® | 
u — ut il). 
(4+»)n:n®+ erNz 


Aus Gl. (2) folgt die durch das Biegungsmoment 
verursachte Querverschiebung 


A n?n? + 8 — 


cos 


\ 


z (Ar) p Wr (+ Felt). 


(Arm mE (24) 


Da die Summe der Teilverschiebungen die ganze Quer- 
verfschiebung ergeben muß, folgt /,=0 A=--f.. 
Demnach ist: die durch das Biegungsmoment verur- 
sachte Querverschiebung » 


2 
nn + 8, 
w=(2-+») 


m w—fı(). 
ern tsßtrN)n 


M= nem 


Soll nicht nur die ganze Querverschiebung, sondern 

'sollen auch die Teilverschiebungen an den Enden 

=( sein, so folgt f, = 0. Kann (4 + v) n?n? gegen 
A 


8(2 +) vernachlässigt werden, also bei den nie- 


deren Frequenzen langer dünner Federn, ist 


I+r..E 
oa=n?!n:D. I@Hr)mB; ET 


l+»v FD: x 


DOT et E 
Ga FE Kucca sin (wt +0) £ 


die durch die Querkraft verursachte Qyerverschiebung 
2 
kann vernachlässigt werden, Kann ze gegen n?n? 


vernachlässigt werden, also bei den höheren Frequenzen 
kurzer dicker Federn, so ist 


ea Vz nr 1l-+» a D* 


Be 7 a er 
4149, NER. 
= Er, FF Ky cos —— sin(wi ta) , 
Ba; Bene a 
a a a 


‚die durch die Querkraft verursachte Querverschiebung 


ist nur wenig kleiner als die durch das Biegungs- 
moment verursachte Querverschiebung. 

6. Vernachlässigung des Summanden Träg- 
heitsmoment mal Winkelbeschleunigung. Für 
die Querverschiebung gilt die Differentialgleichung 

dw m. a A u A ER 
dx: 2(1 +») FI 9220: l+, FID 9: 
Für die Eigenfrequenzen gilt 
4 7g4 
RE a u 2 
nm’ +42 49) 
2 72 
we 2n?n®(1l-+v) Flo; 
n?n?®+4(2 +9) m 


\ r 
g9= un en cos sin (ot+a); 
"art rn; 
Mara. e 
Wy, 1 w; 
n?n? + 4(2 +») Di 
12 
at Nm 
We W. 


n2n2 + 4(2 +9 


ee f 
Kann 4(2--») 7; gegen n?n? vernachlässigt werden, 


so ist 


o=nn y>c« 4 M=—2(l+»)Fiw; 


Q=2nn(l +) Rz 008 
die durch das Biegungsmoment verursachte Quer- 
verschiebung kann vernachlässigt werden. 

7. Vernachlässigung der durch die Querkräft 
verursachten Querverschiebung. Für die Quer- 
verschiebung gilt die Differentialgleichung 

1l+» D* dw m? Aw. '0:w 
tv 20m, 8 dan ae 
Für die Eigenfrequenzen gilt 
I+v, mt, F 


en 
Zee, 


sin (w2 + a); 


=. 


= 4 
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‚Nach Gl. () ist das Bingungsmohienk, 
BREI DE RR re HERZ 
tere el 24» a0 
83 
Q=4 Er ER" wit). 
hi n2. Pen h } l i 


3 
Kann 8 sh gegen n’n? vernachlässigt werden, so ist 


B 1,7 
ernallir. mi’ 


Ih» 
Q gar "FK,0 cos = ® sin (ot +a). 
‚8. Vernachlässigung des Summanden Träg- 
heitsmoment mal Winkelbeschleunigung und 
der durch die Querkraft verursachten Quer- 
verschiebung. ‚Für die Querverschiebung gilt die 

Differentialgleichung 


147,2 tw 0: 


PETE Fra re 
Die Eigenfrequenzen sind 
‘ l1+v F 
mar? N 
o=n 2 nn Im‘ 
ee — 
= —_n?n N 
n ar ——W; 
1-+» FD: 
Ren 
Qa=n mt ppf:i° cos" TE sin(ot + a). 
Hannover. RR Ludwig.. 


Eine Beziehung zwischen der nichtlinearen 
Elastizitätstheorie und‘ der Verfestigungs- 


‚theorie von Ros-Eichinger-Schmidt. 


M.Rosund A. Eichinger[1]!) haben im Jahre 1926 
durch Versuche nachgewiesen, daß im Verfestigungs- 
gebiet eine funktionale Beziehung zwischen ger Inten- 
sität des Spannungsdeviatore und der Intensität des 
Deformationsdeviators folgender Art 


Yo, — om)? + (0, — 0m)? + (0, — om)? \ nr (1) 
= (Ve — em)? + (a — Em)? + (s—Em)?) | 


‚ existiert, die vom Beanspruchungsweg unabhängig ist 


und nur vom Material selbst abhängt. Dieser Tat- 
bestand, der von R. Schmidt [2] durch weitere Ver- 
suche bestätigt wurde, hat: zwar zu einigen Arbeiten 
Anlaß gegeben [3], jedoch ist bis jetzt nicht bemerkt 
worden, daß man daraus ein bemerkenswertes Ergebnis 
ableiten kann,nämlich, daß zwischen der Verfestigungs- 


‘theorie, die auf Gl. d) gestützt ist, und der nicht- 


linearen Elastizitätstheorie kein Unterschied vorhan- 
den ist. Dieses Ergebnis wollen wir hier ableiten. Wir 


'beabsichtigen, in einer späteren Arbeit Anwendungen 


zu geben. 
Wir machen folgende Annahmen: 
1. Der Körper bleibt während der LEER 


'isotrop. Daraus folgt, daß die Hauptdehnungs- 


richtungen dauernd mit den Hauptspannungsrichtungen 
zusammenfallen. 

2. Die Figur der drei Mohrschen Hauptdehnungs- 
kreise in der &—-y-Ebene ist stets ähnlich zur Figur 
der drei Mohrschen Hauptspannungskreise in der 
o—t-Ebene. Dies kann man: so ausdrücken 


EN RE 


Ir Yarshs 
Vorteilhafter für uns-ist es, statt‘ Gl. (2) 
RR Ver, VRR Z Ta LO TRRRBRIT €. 
€&1 Em € Em €; —Em 


zu schreiben. Die zwei Formen (2) und (3) sind mit- 
einander identisch, und gerade das Ersetzen von Gl. (2) 


1) Die Zahlen in eckigen Klammern verweisen auf das 
Schrifttumverzeichnis am Schluß der Arbeit. 


-und mit Hilfe von (5) bekommen wir 
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durch G1. (3) bildet einen wichtigen Schritt für unsere 
Ableitung. 

Die Keckische Elastizitätstheorie werden wir in der 
Weise gewinnen, daß wir außer diesen zwei Blue 
noch folgendes voraussetzen. 

3. Es ist immer 


m a 92 = 0m 4, 0, — Om. 


—Em © h gem - Em & Em 
8 Es ist immer : 


= — =3k = konst, ne dire (0) 


Diese zwei Stoffkonstanten @undksind der Schub- 
modul und der Kompressionsmodul. 
Aus (4) folgt nun 


—m= 26 (e1— Em) 


0, = 26 (61 — em) + 0m = 20, + m (3220| 


0, = = 20|., + 2m (8:20) (6) 
N, — 20 + 37020] 


Diese Gl. (6) sind jene der klassischen Elastizitäts- 
theorie. Die natürlichste nichtlineare Elastizitäts- 
theorie wird man daraus gewinnen, daß man die 
Konstanten @ und k durch Funktionen ersetzt. Man 
wird also 


i-m-hli—em) Ü=1239...0.. Mm 
und ; 
Oom=P% (&m) a a I Er (8) 


voraussetzen, wobei f; und vier Funktionen sind, 
die nur vom Material abhängen. Die Gleichungen der 
‚nichtlinearen Elastizitätstheorie werden also folgende 
"sein: 


°, = fı(8ı— Em) +9 (em) 
0, = fg (eg — Em) + P-(em) En), 


= [3 (es — Em) + 9 (em) 
Wir ee nun, daß aus Gl. (3) Er, folgt: 


—— (10) 
en: (,— Im)? + (09 — Om?) + (03 — om)? 
Y (&ı — Em)” Aus (eg — Em)’ + (&£3 — &m)* , 
Die nichtlinearen  Elastizitätstheorien 
können wir also auch so gewinnen, daß wir 
anstatt 
(ei — om) = Fi (ei — Em) 
zu schreiben, lieber 
V (0, — om)? + (0, — Om)? + (03 — Om)? 
= F (Yleı — em)? + (eg — em)? + (&5— Em?) 
schreiben, was freilich genau dasselbe ist. 
Der Vergleich mit (1) zeigt uns, daß gerade der 


Umstand,daß im Se Pd de die zwei 
Größen - 


(=1,2,3) 


V (&1 — Em)® ze (&9 — Em)? y+ (&s — Em)” 


und 


Y(, — Om)” + (0, — Om)? + (03 — om)? 
undkeineanderen miteinanderin funktionaler 
Abhängigkeit stehen, uns zum logischen 
Schluß bringt, daß zwischen der nichtlinearen 
Elastizitätstheorie und der Verfestigungs- 
theorie, dieauf (l) aufgebaut wird, kein Unter- 
schied existiert. 

Damit ist unsere Behauptung schon bewiesen. Dar- 
über hinaus bemerken wir, daß die Gl. (9) oder auch 
die damit gleichbedeutenden Gleichungen 


9 x Nachrichten 


_ Ve en + Fe]. 
V (&ı Em)? + (&8 — Em)? + (es Em)? 
. (Em) + Plem) 
ER I VER Be 
6 N [Yleı — em)? + (es — em)? + (&3 —em)?] 
/ Ve, — em)? + (&3 — Em)’ + (&3 A Em)” “N 
| (9 — Em) + Flem) 
_ FlVle, — em)? + (eg — em)? + (8 — em)®] 
cr EEE BEE TER VE FE De er TONER) 
Ver — em)? + (e,— Em)? + (es — Em)? 
(Es — Em) + P (em) 
sowohl die Gleichungen der nichtlinearen Elasti- 
zitätstheorie: wie auch jene der Verfestigungs- 
theorie von Ros-Eichinger-Schmidt sind. Zu 
bemerkenist aber.noch, daß in der hier gegebenen Form 


diese Verfestigungstheorie viel weittragender darge- 
stellt ist, als dies. bis jetzt der Fall war, da wir jetzt 


0] 


05 
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auch eine eventuelle 
war, daß die Volumenänderung gleich Null ist. 


Schriftt um, 


1]M. Ros und A. Eichinger: Versuche zur Klärung 
ve der Frage der Bruchgefahr. Eidg. Materialprüfungs- 
anstalt. Zürich 1926 und 19299, 


[2] R. Schmidt: Über den Zusammenhang von 


Spannungen und Formänderungen im Verfesti- 
ö Mean. © Ing.-Arch. Bd.3 (1932), 8. 215 bis 235.. 


3] R. Moufang: Das plastische Verhalten von dünn- 
en wandigen Holinen unter statischem Innendruck. 


Z. angew. Math. Mech. Bd. 20 (1940), & 24. 
München). A. Philippidis. 
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®) Lehrstuhl für technische Mechanik, Technische Hoch- 
schule München. Herrn Prof. Dr. L. Föppl und der 
Förderergemeinschaft der Deutschen Industrie möchte ich 
hier herzlichst danken für die Unterstützung, die sie mir 
bei meiner Arbeit gewährt, haben, 


NACHRICHTEN 


5 Tagung in St. Louis. 

In der Zeit vom 23, bis zum 25. Oktober fand in 
St, Louis eine Tagung statt, auf der Probleme der 
Strömungsforschung behandelt wurden, Es wurden 
dort die folgenden Vorträge gehalten: 


SCHARDIN (Weil a.Rh.) — Bedeutung der Gas- 


dynamik für die Ballistik. E 

SAUER (Weil a.Rh.) — Methoden für die praktische 
Berechnung von Überschallströmungen. 

Dr. HEINZ (Lörrach/Baden) — Fehlerbetrachtungen 
zum Charakteristikenverfahren. 

SCHUBERT (früher Peenemünde) — Theoretischeund 
experimentelle Untersuchungen am axial ange- 
strömten Kreiskegel. 

SAUER (Weila. Rh.) — Die Berechnung der Strömung 
um einen schief angeblasenen Kreiskegel (Diskus- 
sionsbeitrag). ” 

BETZ (Göttingen) — Uber den abgelösten Verdich- 
tungsstoß. N 

WALZ (Göttingen) — Über die Grenzschicht bei hohen 
Machschen Zahlen. 

MANGLER (Göttingen) — Kompressible Grenz- 
schichten ‚bei drehsymmetrischer und bei ebener 
Strömung. 

WENDT (Braunschweig) — Die kompressible laminare 
Grenzschicht an der Platte und am Kreiskegel. 

KETTEL — Turbulenzmessungen mit Hilfe der Ionen- 
sonde, 

-KRAUS — Oberflächentemperatur an schnellfliegen. 

den Körpern. 


REUTTER (Karlsruhe) — Ebene Potentialströmung 
in der Umgebung der Schallgeschwindigkeit, % 
HAACK (Bad Gandersheim) — Theoretische und experi- 

mentelle Untersuchungen der Stabilität. 

ATHEN (Elmshorn) — Übereinige Fragen des Gesamt- _ 
widerstandes von Rotationskörpern in Zusammen-. 
hang mit dem Optimumproblem. : 

KOHLER (Greifswald) — Die gaskinetische Theorie 
des Verdichtungsstoßes. . 


In den Jahren 1945 und 1946 starben von den 
Mitarbeitern der Zeitschrift, für angewandte 
Mathematik und Mechanik: 


Prof. Dr. R. Baldus, München. 

Prof. Dr. W.Cauer, Berlin. 2 
Prof. Dr. H, Föttinger, Berlin. S 
Prof. Dr. R. Fuchs, Berlin. 

Prof. Dr. J. Horn, Darmstadt. 

Prof. Dr. W. v.Koppenfels, Brünn. 
Prof. Dr. G.Kutta, Stuttgart. 

Prof. Dr. M.Lagally, Dresden, 

Prof. Dr, W. Ludwig, Dresden. 

Prof. Dr. OÖ. Mader, München. 

Prof. Dr. L. v. Schrutka, Wien. 

Prof. Dr. Th. Vahlen, Berlin. 
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